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A 12.4 day periodicity in a close binary system after a su- 
pernova 
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Neutron stars and stellar-mass black holes are the remnants of massive star explosions! Most 
massive stars reside in close binary systems”, and the interplay between the companion star 
and the newly formed compact object has been theoretically explored?) but signatures for 
binarity or evidence for the formation of a compact object during a supernova explosion are 
still lacking. Here we report a stripped-envelope supernova, SN 2022jli, which shows 12.4- 
day periodic undulations during the declining light curve. Narrow Ha emission is detected 
in late-time spectra with concordant periodic velocity shifts, likely arising from hydrogen gas 
stripped from a companion and accreted onto the compact remnant. A new Fermi-LAT y- 
ray source is temporally and positionally consistent with SN 2022jli. The observed properties 
of SN 2022jli, including periodic undulations in the optical light curve, coherent Ha emission 
shifting, and evidence for association with a y-ray source, point to the explosion of a massive 
star in a binary system leaving behind a bound compact remnant. Mass accretion from 
the companion star onto the compact object powers the light curve of the supernova and 
generates the 7-ray emission. 


SN 2022jli was discovered by Libert Monard on 2022 May 5 (JD=2459704.67) and was later 
confirmed by several surveys (Method Section § 1). It was classified as a Type Ic supernova (SN 
Ic). The supernova explosion occurred in the spiral arm of a nearby galaxy, NGC 157 (Extended 
Data Fig. (I), at a redshift of z = 0.0055 with a peculiar-velocity-corrected Hubble flow distance 
of D = 22.5Mpc (Method Section § 2). The photospheric spectra of SN 2022jli match well to 
those of spectroscopically regular SNe Id. We measured an ejecta expansion velocity of around 
8,200kms~! from the absorption minima of prominent Fe II lines (Extended Data Fig. [2). Re- 
markably, the supernova brightens again around one month after discovery, which is unusual for 
SNe Ic. Since then, we obtained extensive follow-up photometry and spectroscopy data (Methods 
Section § 3 and § 4). 


The light curve of SN 2022jli (inset panel of Fig.|1) shows three distinct evolutionary phases: 
Phase I, the first decline phase; Phase 1, the rebrightening and gradual decline phase; and Phase III, 


the late-time fast decline phase. During the gradual decline in Phase II starting around two months 
after discovery, SN 2022jli shows periodic undulations in the light curves (Fig. (ip. These bumps 
appear in all the observed optical bands and last around 200 days until the onset of the late-time 
rapid decline. We performed a periodogram analysis of the multiband light curves and found 
a prominent peak at 12.4 days in the power spectrum (Fig. B). The false alarm probability for 
the detected periodicity is less than 107° (Method Section § 5). We also performed a periodicity 
analysis on the individual band light curves and found the same prominent period in each band 
around 12.4 days (Extended Data Fig. B). The phase-folded light curve adopting P = 12.4 days 
(Fig. |2) reveals a constant profile composed of a fast rise lasting ~ 3 days (~ 0.25 phase) and a 
relatively slow decline. We divided the nearly 200-day undulating light curve into two halves with 
equal time spans and repeated the periodicity analysis performed for the whole dataset. We did not 
find significant differences in either the bump period or profile between the two parts. Individual 
bumps have been observed in hydrogen-poor superluminous supernovae” and evidence for rapid 
variability in the optical light curve has also been found in Type Ia SN 2014¥% but this is the first 
time multiple bumps with such a strong periodic signal have been detected. 


Figure [3| shows a selected sample of spectra between +139 and +280 days after discovery. 
Prior to the rapid decline phase (K 270d), these spectra show multiple emission features as well 
as noticeable continuum emission. Compared with the nebular spectra of other stripped-envelope 
supernovae (SESNe), our spectra show stronger iron-group-element emission and many oxygen 
emission lines (Extended Data Figs. M). The comparison reveals that SN 2022jli resembles 
some hydrogen-poor superluminous supernova H] transients that might have an association with 
Gamma-ray Burst (GRB)3) and other peculiar long-lasting supernovae, such as SN 2012au"4! 
and iPTF15dtg!5 which are suggested to be powered by central engines. During the gradual de- 
cline phase, SN 2022jli shows unique emissions around 6500 A with a relatively narrow feature 
superposed on a broad component. We attribute the narrow feature to hydrogen Ha emission 
(Method Section § 9). The Ha lines show remarkable shifting around zero velocity in a repetitive 
pattern consistent with having a 12.4-day period (bottom panels of Fig.[3p. SN 2022jli experiences 
significant spectroscopic evolution as it enters the fast decline phase (Extended Data Fig. [5). One 
conspicuous change is that the relatively narrow emissions, for example, the permitted O1 lines 
with a full width at half maximum (FWHM) around 2,600 kms~, disappear, and broader [O 1] 
6300, 6363 with an FWHM around 5500 kms~! emerge (Extended Data Fig. |8). 


The pseudo-bolometric light curve of SN 2022jli (Method Section § 6) displays an evolution 
similar to those of the individual bands. It displays two peaks (Fig. pp. making it a double-peaked 
SESN, though the rising part of the first peak is missing. The Phase I light curve could be explained 
by °°Ni radioactive decay as in normal SNe Ic (Method Section § 7; see discussion below). The 


Phase 11 light curve requires another energy source of around 2 x 10“ ergs, and the dramatic drop in 
Phase III implies a sudden shutoff of the extra energy input. To probe the properties of such a late- 
time energy source, we obtained two epochs of radio observation with the Australian Telescope 
Compact Array (ATCA) at +213 and +228 days after the discovery, respectively. Neither epoch 
had a clear detection of SN 2022jli, and the second epoch with the longer exposure time gave 
5c limits of < 0.074mJybeam™! and < 0.055mJy beam! at central wavelengths of 5.5 GHz 
and 9.0 GHz, respectively. We also performed X-ray observations of SN 2022jli with NuSTAR 
(+227 to +237 days after discovery) and the Chandra X-ray Observatory (+257 to +266 days 
after discovery), but did not detect any emission in soft or hard X-rays (Methods Section § 11) 


with upper limits of L39_go Kev < 2.5 x 10*° ergs! and Lo.5_7Kev < 1.3 x 10° ergs}. 


We searched for high-energy y-ray emission using data from Large Area Telescope on board 
the Fermi Gamma-ray Space Telescope (Fermi-LAT ) and found a y-ray source in the direction of 
SN 2022jli (Method Section § 15). The new y-ray source is detected after the supernova explosion, 
and there was no detection in the past 13.5 years of archival data before the supernova explosion 
(Extended Data Fig. ug SN 2022jli is within the 68% confidence localization area of the new 
source, and the detection time of the y-ray photons from the new source shows evidence for a cor- 
relation with the 12.4-day optical flux undulation (Extended Data Fig. (8). The temporal and spatial 
coincidence, together with potential periodicity, suggests an association of the new y-ray emission 
with SN 2022jli. The new source is most significantly detected in the 1-3 GeV energy band with 
Lisey = 3.1 x 10*! ergs~! in November and December 2022, and the whole y-ray light-curve 
is shown in Fig. |4| We suggest the y-ray emission was only detected several months after the 
supernova explosion due to high pair-production opacities to y-ray photons at early time! The 
observed y-ray emission falls below the Fermi-LAT sensitivity at the end of December 2022, one 
month prior to the fast drop in optical flux. 


Below we discuss the energy source that powers the unique light curves of SN 2022jli. Due 
to the photometric and spectroscopic resemblance to normal SESNe in Phase 1 (Extended Data 
Figs. [2|and|9p, it is natural to attribute the first peak to the same origin as in normal SESNe, a °°Ni- 
decay powered peak. It is more intriguing to think about the energy source of the second peak 
and the powering mechanism of the undulations. Supernova ejecta interaction with hydrogen- 
poor circumstellar medium (CSM) has been commonly adopted to explain the observed bumps 
in the light curves of SESNeSP118120 It is appealing to connect the periodic bumps to ejecta-CSM 
interaction (ECI) with a CSM having density fluctuation, such as the nested dust shells surrounding 
the Wolf-Rayet binary WR 14024) and a late-time drop in the light curve is also expected after the 
ejecta sweeps through the confined dense CSM2223 However, the persistence and short period 
of the bumps, especially the short time scale of around three days for the rising part, put tight 


constraints on the feasibility of the ECI scenario to explain the observed properties. For example, 
consider a CSM composed of nested spherical shells. In that case, the light travel time difference 
á š R Vej t 
between the near and far sides of the ECI is At ~ + = 10 X torkos 300 days 
effective ejecta expansion velocity and t is the time from the explosion. For the observed ejecta 


days, where vej is the 


velocity of 8,200 km st, at 260 days, the light travel time difference is around 7 days, which would 
smear out any periodic signal from ECI due to CSM density fluctuation. Moreover, the light travel 
time difference is phase-dependent as the ejecta expand and interact with CSM at different radii. 
The observed periodic bumps span around 200 days, showing no significant evolution, arguing 
against the ECI explanation for the periodic undulations. 


The observed bumps have a short time scale, a constant profile, and a constant ratio to the 
continuum flux, i.e., the fluctuating light is not being diluted even when the ejecta expands. This 
requires an emitting region that covers a constant fraction of the expanding ejecta and that behaves 
coherently. The corresponding power source must reside in a spatially confined region to remain 
coherent and be located close to or at the center of the ejecta. In particular, the quick rise of the 
bump profile requires that the diffusion time of optical photons through the ejecta is short (< 3d). 


The periodic undulation first appears around the second peak, requiring that the diffusion 
time of the ejecta is already short at this time. This means that the Phase II light curve could not be 
powered by °°Ni decay because the diffusion time would be similar to the rise time to peak (tens of 
days), which is inconsistent with powering the much shorter bumps. This argument applies to any 
other energy source that was produced during the explosion (such as a magnetar) with a decreasing 
or constant energy output after the explosion. A possible mechanism that could provide extra 
energy to power the later evolution of this supernova is accretion onto a newly-formed compact 
object2422) which we propose to explain the observed light curves of SN 2022jli (Method Section 
§ 7). In this model, the deferred onset of the energy input could be due to the time it takes to 
form an accretion disk, while the abrupt luminosity drop-off at the late time is accounted for by the 


central object running out of infalling gas to fuel the accretion. 


Given the periodic signals observed in three different aspects of the supernova, the undulation 
in the optical light curve, the shifting Ha emission, and the high-energy y-ray photons, the most 
viable explanations seem to inevitably involve a binary system to provide the “clocking” mecha- 
nism. The luminosity of the narrow Ha emission closely follows the evolution of the bolometric 
luminosity (Extended Data Fig. (10), even during the transition phase when the supernova lumi- 
nosity shows an enhancement before diving down. This means the Ha emission shares the same 
origin as the excess luminosity of the supernova, implying the Ha also comes from the central 


region of the system. Therefore, a companion star with a hydrogen-rich envelope is necessary to 


provide the hydrogen, which corroborates the existence of a binary system. Such a binary system 
provides direct evidence supporting the binary origins of some SESNe. 


We suggest that the compact supernova remnant and the companion star remain bound in 

a binary system after the supernova explosion. One plausible idea to power the supernova is 

through the accretion from the bloated companion star onto the compact object, forming an ac- 

cretion disk. An eccentric orbit is expected (Method Section § 12), which modulates the accretion 

rate generating the observed undulation. With L = eMe, and adopting the optical-NIR pseudo- 

bolometric luminosity for the accretion luminosity L, we calculated a peak accretion rate M around 
€ 


< =1 A 
4 x 1078 (z) Mz yr™t, and a total mass to account for the accretion-powered energy around 
& 


-1 : MOP : . . 
1073 (z) Mo, where e is the radiative efficiency of accretion. Such an extreme accretion rate 
could give rise to an extreme outflow of wind and/or a jet. The generation of a jet and/or wind 
favors the idea that the detected y-ray emission of SN 2022jli shares a similar origin to that seen 


in black hole X-ray binary systems“ Jets have been predicted in several supernova explosion 
scenarios2132| 


The unprecedented properties of SN 2022jli tell that whatever happens in the system should 
be a rare phenomenon, which might be explained by the rarity of a bound binary system surviving 
a supernova explosion B3B4 SN 2022jli provides direct observational evidence for the survival 
of such a binary system after a supernova explosion. SN 2022jli builds a direct link between 
the supernova explosion and the formation of a compact object. It highlights the importance of 
the interplay between the companion star and the newly formed compact object in shaping the 


appearance of a core-collapse supernova. 
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Figure |1| Multiband light curves of SN 2022jli showing periodic undulations. (a) g,7,7, z,c, 
and o-band light curves of SN 2022jli during phases between +50 and +260 days after discovery. 
The dashed lines show the polynomial fit to the light curve, which serves as the “baseline”. An 
empirical model with a 12.4-day period is shown as a solid line for each band. The inset panel 
shows the whole range of the multiband light curves. (b) The relative undulations in g,1r, 7, z,c, 
and o-band light curves. All the error bars of the data points are 1o confidence intervals. 


Figure|2} Multiband periodogram and the undulation profile of SN 2022jli. (a) Power spectrum 
of the multiband light curve by jointly fitting of g,7,2,z,c, and o-band light curves shown in 
Fig.|1| The power spectrum shows clear peaks at 12.4 days and lower-frequency harmonic aliases 
at around 24.8, 37.2, and 49.6 days. (b) The folded undulation profile using the period of 12.4 days. 
The black dashed lines show the best-fit empirical model (Section § 5) describing the undulation 
profile by a fast rise and gradual decline. Error bars of the data points are 1o confidence intervals. 


Figure |3} The spectral evolution of SN 2022jli between +139 and +280 days after discovery 
(a) Flux calibrated optical spectra of SN 2022jli. The spectral phase relative to the supernova 
discovery time is shown with the color bar. (b) The evolution of Ha emission. The left panel shows 
the zoom-in view around 6000-6800 A. The right panel shows the narrow Ha after subtracting the 
pseudo-continuum (indicated by the dash-dotted line in the left panel). For each spectrum, the 
observation date is shown on the left, and the phase relative to the peak of the bump profile is 
shown on the right. The Ha velocity is shown at the top. (c) The same as the right panel of (b), but 
the flux of each spectrum has been scaled to have the same integrated flux within 6430—6680 A, 
and the spectra are sorted by the phase relative to the minimum of the bump profile. 


Figure [4} The pseudo bolometric light curve and multi-frequency data of SN 2022jli. (a) The 
y-ray source detection map generated with bimonthly Fermi-LAT observation in the energy band 
of 1—3 GeV using the Poisson noise matched filter method (Section § 15). The red plus symbol in 
the center of the field indicates the position of SN 2022jli. A clear y-ray source is detected in the 
2022 November 1 to 2023 January 1 bin. (b) The multi-frequency light curve of SN 2022jli. The 
black points show the pseudo-bolometric light curve from 3750 A to 25000 A. The blue line shows 
the radioactive decay model with 0.15 Mọ °°Ni. The light curve of high-energy (1 —3 GeV) y-ray 
emission associated with SN 2022jli is shown with blue points. Two epochs of X-ray observations 
with NuSTAR (yellow) and Chandra (magenta) around 250 days resulted in non-detections. All 
the error bars are lo uncertainties. All the non-detections are shown as 30 upper limits. 
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Figure 1: Multiband light curves of SN 2022jli showing periodic undulations. (a) g,r, 7, z,c, 
and o-band light curves of SN 2022jli during phases between +50 and +260 days after discovery. 
The dashed lines show the polynomial fit to the light curve, which serves as the “baseline”. An 
empirical model with a 12.4-day period is shown as a solid line for each band. The inset panel 
shows the whole range of the multiband light curves. (b) The relative undulations in g,7, 72, z,c, 
and o-band light curves. All the error bars of the data points are 1o confidence intervals. 
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Figure 2: Multiband periodogram and the undulation profile of SN 2022jli. (a) Power spectrum 
of the multiband light curve by jointly fitting of g, r, i, z, c, and o-band light curves shown in Fig.[1| 
The power spectrum shows clear peaks at 12.4 days and its lower-frequency harmonic aliases at 
around 24.8, 37.2, and 49.6 days. (b) The folded undulation profile using the period of 12.4 days. 
The black dashed lines show the best-fit empirical model (Section § 5) describing the undulation 
profile by a fast rise and gradual decline. Error bars of the data points are 1o confidence intervals. 
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Figure 3: The spectral evolution of SN 2022jli between +139 and +280 days after discovery. (a) 
Flux calibrated optical spectra of SN 2022jli. The spectral phase in days relative to the discovery 
time is shown with the color bar. (b) The evolution of Ha emission. The left panel shows the 
zoom-in view around 6000-6800 A. The right panel shows the narrow Ha after subtracting the 
pseudo-continuum (indicated by the dash-dotted line in the left panel). For each spectrum, the 
observation date is shown on the left, and the phase relative to the minimum of the bump profile is 
shown on the right. The Ha velocity is shown at the top. (c) The same as the right panel of (b), but 
the flux of each spectrum has been scaled to have the same integrated flux within 6430-6680 A, 
and the spectra are sorted by the phase relative to the minimum of the bump profile. 
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Figure 4: The pseudo-bolometric light curve and multi-frequency data of SN 2022jli. (a) The 
y-ray source detection map generated with bimonthly Fermi-LAT observation in the energy band 
of 1-3 GeV using the Poisson noise matched filter method (Section § 15). The red plus symbol 
indicates the position of SN 2022jli. A clear y-ray source is detected in the 2022 November 1 to 
2023 January 1 bin. (b) The multi-frequency light curve of SN 2022jli. The black points show the 
pseudo-bolometric light curve from 3750 A to 25000 A. The blue line shows the radioactive decay 
model with 0.15 Mọ °°Ni. The light curve of high-energy (1-3 GeV) y-ray emission associated 
with SN 2022jli is shown with red points. Two epochs of X-ray observations with NuSTAR (yel- 
low) and Chandra (magenta) around 250 days resulted in non-detections. All the error bars are lo 
uncertainties. All the non-detections are shown as 30 upper limits. 
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Methods 


1. Discovery of SN 2022jli: SN 2022jli was first discovered and reported to the Transient Name 
Server (TNS; by Libert Monard’? on 
2022 May 5. It was later recovered by wide-field transient surveys including the Asteroid Terrestrial- 
impact Last Alert System (ATLAS2 ATLAS22oat), the Gaia transient survey! (Gaia22cbu), the 
Panoramic Survey Telescope and Rapid Response System (PanSTARR SE PS22gwo), and the 
Zwicky Transient Facility (ZTE EJ, ZTF22aapubuy). The transient was discovered as it was 
rising in the morning sky. Solar conjunction before the discovery of SN 2022jli hindered any 
direct constraint on the explosion time so we set the discovery epoch as the reference time for 
phase definition whenever the reference epoch is not explicitly specified throughout this paper. 
The Gaia detection gives a transient coordinate of right ascension a = 0034™45.690° and dec- 
lination 6 = —08°23'12.16” (J2000.0). The Galactic extinction in the direction of SN 2022jli is 
E(B — V)mw = 0.039 mag. The latest non-detection was on 2022 February 6, from Gaia, 87.5 
days before the discovery. Therefore, due to the large seasonal gap, we barely have any constraint 
on the explosion time of SN 2022)li. 


2. Host galaxy and extinction: SN 2022jli is located on a spiral arm of NGC 157 (Extended Data 
Fig. (ip. NGC 157 has a large number of HII regions, resulting in a complex background with 
strong nebular emission from the host galaxy, which commonly causes over-subtraction of host- 
galaxy lines. The distance to NGC 157 is uncertain. The peculiar-velocity-corrected Hubble flow 
distance gives D (Virgo + GA + Shapley) = 22.5 + 1.6 Mpc (u = 31.76 mag)! using a standard 
ACDM cosmology with Qm = 0.27, Qa = 0.73 and Ho = 73 km s~! Mpc7! ®©! The Tully- 
Fisher method gives smaller values of 12-13 Mpita In this work, we adopt the distance of 22.5 
Mpc for our analysis but caution that distance-related quantities are subject to large uncertainties. 


We detected prominent sodium absorption from the host galaxy in the supernova spectra, 
distinguished as a doublet in our medium resolution spectra taken with IMACS and X-Shooter 
(Extended Data Fig. [I). We measured a sodium absorption equivalent width EW(NaI D1+D2) = 
1.07 + 0.11 A, which converts to E(B — V)nost = 0-25? 4 mag E We did not detect significant 
sodium absorption from the Milky Way, consistent with the low Galactic reddening. We also 
detected narrow diffusion interstellar band (DIB) absorption at 6283 A (Extended Data Fig. lib, 
from which we measured an equivalent width of EW(DIB6283) = 0.27 +0.03 Å, which converts to 
E(B —V rost = 0.25+40.03 mag“®l or E(B—V aost = 0.28 40.04 mag“?! The extinction uncertainty 
derived from DIB absorption only includes the uncertainty of the equivalent width measurement 
but not the uncertainty from the conversion. In this work, we adopted E(B — V nost = 0.25 mag 
and the conservative uncertainty obtained from sodium absorption for the analysis. 
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3. Photometry: SN 2022jli was detected by ZTF starting from 2022 June 27. ZTF gri photometry 
obtained with the ZTF survey = camera was processed with the ZTF image reduction pipelind=! 
employing the ZOGY image-subtraction method} We obtained additional gri images with the 
robotic 60-inch telescope at Palomar (P6033), using the Spectral Energy Distribution Machine 
(SEDM*4). FP i pel was used to extract PSF photometry from image subtraction against Sloan 
Digital Sky Survey (SDSS) templates. ZTF and SEDM photometry were obtained and reduced in 
real-time mode and streamed to the Fritz SkyPortal 5557 to aid further follow-up observations. The 
SEDM photometry obtained above with FPipe was considered preliminary and was refined as 
described below. 


We included photometry data reported to TNS by Libert Monard at the Kleinkaroo Obser- 
vatory (KKO), which were taken on 2022 May 5, 6, 11, and 22. We obtained All-Sky Automated 
Survey for Supernovae (ASAS-SN) g-band photometry * from the ASAS-SN Sky Patrol 59] We 
adopted the option of “Image Subtraction (No reference flux added)” for the photometry method, 
which performs aperture photometry on the coadded image subtracted data for each epoch but does 
not add the flux of the source on the reference image to the light curve. We noted a constant offset 
between ASAS-SN g-band and ZTF g-band photometry in flux. We added a constant flux of 360 
mJy to the measured flux in the ASAS-SN data to minimize the difference between ASAS-SN 
and ZTF in the overlapping time range. We extracted ATLAS-c and ATLAS-o band light curve 
data from ATLAS forced photometry server Bee We also used three epochs of Gaia-G band 
photometry obtained by Gaia transient survey E7 through the Gaia Alert service ©) 


After discovering the periodic bumps, we started more follow-up observations with the 0.8m 
RC32 telescope operated by the Post Observatory (PO). We performed aperture photometry on the 
reference-subtracted images and calibrated them with the SDSS standard catalog. For image sub- 
traction, we built template images from SDSS images The FP ipe pipeline did not successfully 
reduce all P60 SEDM images. We performed the same photometry for SEDM images with some 
manual assistance. The detailed photometry procedures are described in ref4l 


We also observed SN 2022jli with the Multi-channel Photometric Survey Telescope (Mephistd®15) 
during the telescope commissioning phase from 2022 December 10 to 2023 February 6. The ob- 
servations were done in wv and iz bands. We conducted image subtraction and then performed PSF 
photometry on the subtracted images. The images taken on the night of 2023 February 6 were used 
as template images. We obtained the flux of SN 2022jli with PSF photometry from the template 
images and added to the flux obtained from the difference images. 


Supplementary Information Fig. |1| shows the observed light curves, and all photometry is 
listed in Supplementary Information Table|I] 
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4. Spectroscopy: We obtained 46 low-to-medium resolution (R ~ 100 to 6000) spectra of SN 2022jli, 
taken with the instruments listed in Supplementary Information Table [2] The spectra are shown in 
Supplementary Information Figs. and f] All spectra will be made publicly available through 
the Weizmann Interactive Supernova Data Repository (WISeREP®I&), Detailed information on 
the spectroscopic observations and data reduction is listed below. 


P60/SEDM: The Spectral Energy Distribution Machine (SEDM is an integral field unit 
(IFU) spectrograph mounted on the 60-inch robotic telescope (P60 =) at Palomar Observatory. We 
conducted 25 epochs of spectroscopy observation with the SEDM between 2022 June 29 and 2023 
February 2. The SED Machine has a very low resolution (R ~ 100) covering the wavelength range 
from 3650 to 10000 A. All SEDM IFU data were reduced using the pipeline described in 70) and 
new modules for the SEDMachine described in were used to remove contamination from cosmic 


rays and non-target light. 


NOT/ALFOSC: The Alhambra Faint Object Spectrograph and Camera (ALFOSC) is mounted 
on the 2.56m Nordic Optical Telescope (NOT). We acquired six epochs of low-resolution spectra 
with NOT/ALFOSC between 2022 August 6 and 2023 February 8. The spectra were obtained 
with a slit width of either 10 or 1’3 depending on the seeing, and using grism #4. The data 
were reduced using the pipeline f oscguiľ The reduction includes cosmic-ray rejection, bias 
corrections, flat fielding, and wavelength calibration using HeNe arc lamps imaged immediately 
after the target. The relative flux calibration was done with spectrophotometric standard stars 
observed on the same night or nights before the observation. 


P200/DBSP: The Double Beam SPectrograph (DBSE®) is mounted on the 200-inch Hale tele- 
scope at Palomar Observatory (P200). The DBSP uses a dichroic (at 5500 A for the used D55 
dichroic) to split light into separate red and blue channels (“sides”), observed simultaneously. We 
obtained three epochs of spectroscopy of SN 2022jli with DBSP. The observations were taken us- 
ing a blue grating with 600 lines per mm blazed at 4000 A, a red grating with 316 lines per mm 
blazed at 7500 A, and a 1”5 wide slit on 2022 August 20 and 2022 November 2, and a 10 wide 
slit on 11 2022 October 11. The data are reduced using the python package DBSP_DRP“ that is 
primarily based on Pype It 214) 


MMT/BINOSPEC: We obtained 5 spectra of SN 2022jli between 2022 December 27 and 2023 
January 31 with Binosped on the MMT Observatory 6.5m telescope. All data were acquired with 
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a grating of 270 lines/mm and a 10 slit mask. The Binospec spectra have a wavelength coverage 
of 3900 — 9240 A. The basic data processing (bias subtraction, flat fielding) is done using the 
Binospec pipeline The processed images are downloaded from the MMTO queue observation 


data archive. All the spectra are reduced with IRAF, including cosmic-ray removal, wavelength 
calibration (using arc lamp frames taken immediately after the target observation), and relative flux 
calibration with archived spectroscopic standards observation. 


Magellan/IMACS & FIRE We observed SN 2022)li with the Folded-port InfraRed Echellette 
(FIRE spectrograph on 2022 August 25, and with both the FIRE and the Inamori-Magellan 
Areal Camera and Spectrograph (IMACS®%) on 2022 December 14. Both FIRE and IMACS are 
mounted on the 6.5m Magellan Baade telescope. The first epoch of the FIRE spectrum was taken 
with the long-slit mode and the second with the echelle mode. The long-slit mode FIRE spectrum 
was reduced with the IDL pipeline £f irehose™| The echelle FIRE spectrum was reduced with 
Pypelt 75176, The IMACS spectrum was taken with the 1200 lines/mm grating at two different 
tilt angles, covering 5130 — 6780A and 7290 — 8920 A. The IMACS spectrum was reduced with 
IRAF in the same way as for the BINOSPEC spectra. 


VLT/XSHOOTER: We obtained three intermediate resolution spectra with the X-shooter echelle 
spectrograph B! on 2023 January 14, 18, 21 through a DDT program (Program ID: 110.25A6, P.I.: 
P. Chen). These data were executed in ToO mode in order to spectroscopically monitor the super- 
nova during the different phases of the light curve undulations. All observations were performed 
in nodding mode and with 1/3/1"2/1"2 wide slits (UVB/VIS/NIR). The observations covered the 
entire spectral range of the X-shooter spectrograph from 3000 to 24800 A. We first removed cos- 
mic rays with the tool astroscrappy4) which is based on the cosmic-ray removal algorithm 
by B3 Then the data were processed with the X-shooter pipeline v3.3.5, and the ESO workflow 
engine ESOReflex #84 A11 data from three arms were reduced in nodding mode. The nodding 
mode reduction is critical for NIR data to ensure a good sky-line subtraction. The spectra of the 
individual arms were stitched by averaging the overlap regions. The atmospheric absorption in the 
VIS and NIR arms was corrected with the software tool mole fit® (v4.2.3). 


5. Periodic undulation in the light curve 


Multiband periodicity analysis We adopt the multiband periodogram method B3, a general ex- 
tension of the well-known Lomb-Scargle approach®221) to quantitatively detect the periodic signal 
in the multiband light curves of SN 2022jli. The light curves in each band are modeled as arbi- 
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trary truncated Fourier series with the period and/or phase shared across all bands. The model 
contains two parts: 1) an Npase-term truncated Fourier series that models the shared variability 
among all six bands, i.e., the “base model”; 2) a set of Npana-term truncated Fourier series for an 
individual band, which models the residual from the base model. The model can be described as 
follows yx(t|w, 0) = A + ae [02 —-1Sin(nwt) + 2,cos(nwt)| + oP + D 0% sin(nwt) + 
oP cos(nwt)] where 6p is the constant offset of the base model, oP is the residual component 
of each filter k, and [02n-1 (0r oP), O2n(Or ot) = [A cosġ, A sing] describes the amplitude and 
phase of the single-component sinusoidal model of d(t) = A sin(wt + ¢). We used the gat spy” 
tool to fit the g, r, i, z, c, and o band light curve observed between +50 and +250 days after dis- 
covery. In practice, the light curves have been detrended and normalized in flux before periodicity 
analysis. We first subtracted fluxes in each band, F, with the mean flux light curve and then 
normalized the residual flux by the mean flux. The mean flux in each band was modeled with a 
polynomial “baseline”, Fy base- The detrending and normalization procedure can be described as 
w = (Fx — Fy base) /F),base- The adopted polynomial orders for the baseline model in each band 
depend on the baseline’s length and smoothness. As a result, polynomial orders of 6, 6, 4, 4, 6, and 
5 have been used for g, r, 7, z, c, and o bands, respectively. We adopted Npase = 4 and Npana = 2 
in our fitting and obtained a significant peak in the power spectrum corresponding to 12.4 days 
which is accompanied by three low-frequency harmonic aliases. We applied the bootstrap method 
to study the significance of the detected periodicity by shuffling the data points and calculating the 
power spectrum as done for the original data. We performed 10° experiments, and none of them 
revealed any similar peaks in the power spectrum, which means the false alarm probability (FAP) 
of the 12.4-day period is smaller than 107°. To test whether there is a significant phase-dependent 
periodicity evolution, we divided the above data into two parts to perform the same analysis on 
each part as done for the whole data. We obtained a best-fit period of 12.23 + 0.10 days for the 
first-half data and a best-fit period of 12.40 + 0.08 days for the second-half data. 3 o uncertainties 
are reported. No significant period change was detected in our data. 


Empirical model of the undulation profile The relative undulation of SN 2022jli (Fig. [2) shows 
characteristic profiles composing a fast rise and then a gradual decline. We constructed an empir- 
ical model for the undulation profile to more quantitatively characterize its features. The basis of 


the empirical model is a piecewise function with two linear components, as given below: 


ce ae 0< t< tise 
w=] ae — (1) 


=A (t-P)+C, tis <t<P 


P—trise 
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where trise is the time for the rising phase, A is the whole amplitude from the lowest to the highest 
point, P is the period, and C is a constant. The empirical model is defined within the range of one 
period t € [0, P], which starts from the point with the lowest value of 4%. To smooth the model, 
we convolved the above function with a narrow Gaussian kernel (o = 0.2 days). We applied the 
above empirical model to fit the phase-folded undulation profile in each band and obtained a rising 
phase duration of 3.5 + 0.2, 2.5 + 0.2, 2.5 + 0.2, 2.7 + 0.3, 1.9+0.2, and 1.9 + 0.2 days for g, r, 
i, Z, €, and o band, respectively. The best-fit empirical models are shown with black dashed lines 
in Fig. [2]and are also used in Fig. [I] 


Individual band periodicity analysis We also performed periodicity analysis with the Lomb- 
Scargle method®224 in individual bands. The left panels of Extended Data Fig. [3]show the Lomb- 
Scargle power spectrum for each g, r, i, c, o, and z band with blue lines. We also compute the 
power spectrum of the window function shown with the yellow lines. The diurnal peaks of the 
window power are pronounced, which is typical for ground-based data with nightly observations. 
The strong diurnal component from the window function causes each frequency signature fo to be 
partially aliased at fọ + nô f, for integers n and ôf = 1 cycle day~'. A significant peak of around 
12.4 days is detected in each band, as shown in the inset panel, where the period corresponding 
to the peak power is given P, = 12.5 days, P, = 12.4 days, P; = 12.4 days, P. = 12.8 days, 
P, = 12.2 days, and P, = 12.2 days. All the detected peaks have FAP smaller than 107°. The 
power spectra of the observational window functions do not show any significant peaks around the 
detected periods, confirming the periodicity’s authenticity. The right panels show the phase-folded 


light curves after detrending and normalization. 


6. Bolometric light curve: We built the pseudo-bolometric light curve with the photometry data 
(Method § 4) and spectral sequence (Methods § 4) of SN 2022jli. We first get the integrated flux in 
the wavelength range of 3750 — 9150 A, which most of our optical spectra cover. This wavelength 
range corresponds roughly to the BV RI bands, and we denote the luminosity obtained from such 
a wavelength range as Lpy rz. Then we derive the fraction of NIR contribution (9150-25000 A) to 
the bolometric luminosity which is defined as fyi = — and apply the correction to get 
Loo = Levrit+Llnir. The U-band contribution to the bolometric luminosity is not considered, but 
is known to be small for SESNe, around 10% around peak light and then quickly declining to less 
than 5% "3! To get the integrated optical flux, for each epoch of photometry, the closest spectrum 
was used as the spectral energy distribution template and scaled by a constant value to match the 
observed flux, then the scaled spectrum was corrected for the Milky Way foreground extinction 
with E(B — V)mw = 0.039 mag" and host extinction with E(B — V )host = 0.25 mag (Methods 
§ 2) with the CCM extinction law adopting Ry = 3.1%! For the NIR fractional contribution, we 
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got fiora = 0.28 and fr204q = 0.29 for SN 2022jli. We derived the same ratios for two epochs 
of spectra of the SESN SN 2013ge and found fsn2013¢e,+94 = 0.24 and fN2013ge,+204a = 0.36. In 
the end, we used a uniform faction of 0.3 for the NIR contribution and added a 10% uncertainty 
to account for the variation. The distance uncertainty was not included in the uncertainty budget. 
The resulting bolometric light curve is shown in Fig. /4]and Extended Data Fig. [6] 


The pseudo-bolometric light curve of SN 2022jli shows clear undulations as seen in the in- 
dividual optical light curves. We constructed the empirical undulation profile the same way as 
done in Method Section § 5 by adopting a 6th-order polynomial for the baseline. We notice an 
enhancement in the pseudo-bolometric luminosity during the last bump prior to the rapid decline 
(panels b and c of Extended Data Fig. (6). If we assume the last bump followed the evolution as 
the previous undulations, the integrated luminosity of the extrapolation in the last bump would be 
~ 3.0 x 10% ergs while the integrated luminosity of observed data in the last bump is ~ 3.9 x 1047 
ergs. In the accretion-powered scenario (Method Section § 7), the excess radiation in the last bump 


€ 


0.01 
of accretion. The undulation profile of the last bump is significantly different from the previous 


j —1 $ F B 
converts to an accretion mass of M excess = 5x 107° ( ) Mo, where e€ is the radiative efficiency 


undulations, with a peak time delayed by around 4 days. 


7. The energy source of the SN luminosity In the standard model of SNe Ic, the SN luminosity 
is dominantly powered by radioactive decay of unstable isotopes, in particular °°Ni, and its daugh- 
ter element °°Co, with some contribution from previously stored kinetic energy in the early time 
through shock cooling emission. During the first month after discovery, SN 2022jli looks like a 
normal Type Ic supernova both in terms of the spectral evolution (Extended Data Fig.|2) and from 
the light curve evolution (Extended Data Fig. (9). Due to the spectral and photometric similarities, 
it is natural to attribute the first decline phase to °°Ni decay as in normal SNe Ic. We found that 
around 0.15 Mọ °°Ni produced in the explosion combined with a characteristic diffusion time scale 
of around 17.5 days and a y-ray escape time scale of around 100 days can explain the first peak 
of SN 2022jli (Fig. (4p. The y-ray escape time scale is adopted arbitrarily with a typical value for 
SNe Ib/c 86. The characteristic diffusion time is a function of ejecta mass, velocity, and opacity. 
Due to the limited amount of data around the first peak, we do not attempt to constrain the ejecta 
properties by modeling the light curve. 


After the first decline phase, SN 2022jli brightens again around +20 days after discovery and 
reaches the second peak around +50 days after discovery. This makes SN 2022jli a double-peaked 
SESN. After the second peak, the luminosity of SN 2022jli declines around ten times from +70 
to +270 days. After subtracting the underlying radioactive decay emission from 0.15 Mọ *°Ni, the 
integrated energy is 2 x 10*° erg, which requires an additional energy source. Double-peaked light 
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curves have been observed in other SESNe, for example, SN 2005bf22, PTF11mnb®! SN 2019cad™ 
SN 2019st LE, SN 2021uy y% SN 2022xx E% and a sample of such objects in ref The mor- 

phology of those double-peaked light curves shows a large diversity. Except SN 202luvy, all of 

them do not have a long-lasting gradual decline phase after the second peak as seen in SN 2022jli 

(Extended Data Fig. |9). The popular explanation for the double-peaked light curves of SESNe in- 

cludes interaction between the ejecta and either extended material at the outskirts of the progenitor 

or detached circumstellar material surrounding the progenitor 1% double-peaked distribution 

of radioactive °°Ni27%, and enhanced magnetar power] Below we discuss the possible origins 

of the second peak of SN 2022jli, and test whether they can explain the undulations observed after 

the second peak. 


The light curve of SN 2022jli after the second peak declines at a rate similar to that of light 
curves powered by fully trapped y-ray from the decay of °°Co, i.e., around 0.01 mag/day. It is 
appealing to attribute the double-peaked light curves to the result of the double-peaked distribution 
of °°Ni. In this scenario, the first peak is powered by °°Ni carried out by a jet-like phenomenon and 
deposited in the outer layers of the ejecta!) and the second peak is powered by °°Ni residing in the 
deep layers of the ejecta. One problem with the °°Ni explanation is the plummet of luminosity at 
very late times. The luminosity declined by one order of magnitude during around 20 days, which 
can not be explained by, for example, y-ray photon leakage or the formation of cold dust in the 
ejecta. Another issue with °°Ni comes when considering the diffusion time scale. The timescale 
of the °°Ni-powered peak is determined by the time it takes the emission to diffuse through the 
ejecta. The long time before the second peak means there must be a large ejecta mass, causing a 
long diffusion time. Such a long diffusion time will smear out short-time scale signals generated 
from the center of the ejecta (Method Section § 12). The other possibility for the energy source is 
the spin-down energy of the newborn magnetar. Similarly, the magnetar model also has problems 
with the late-time drop of the supernova flux and with the diffusion time. 


The interaction between ejecta and CSM could explain the rise and fall of the second peak 
if the CSM is distributed in a confined distance range from the progenitor. The interaction could 
generate periodic energy input if the CSM has evenly distributed density fluctuation. However, 
such periodic energy fails to produce the observed bumps in the light curves of SN 2022jli due to 
the light travel time difference (Method Section § 12). 


Accretion onto the supernova remnant has been proposed as another energy source to power 
supernova light curves, and most of the models consider the fallback ejecta as the fue PE] Ac- 
cretion power could provide the extra energy of SN 2022jli. In light of the existence of a bloated 
companion star in a close orbit (Method Section § 13) and the requirement of the hydrogen-rich 
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material to explain the observed Ha emission (Method Section § 9), we propose that the compan- 
ion star with hydrogen-rich envelope is the donor to fuel the accretion. The accreted mass required 


€ 


0.01 
of accretion. In the accretion scenario, an extremely high Eddington luminosity ratio is inevitable 


to account for the extra 2 x 10*° erg is Mace = 107? ( j Mo where e is the radiative efficiency 
to explain the observed light curve. For a 1.4 Mọ neutron star, the Eddington luminosity ratio 
is around 10*. For a 10 Mọ black hole, the Eddington luminosity ratio decreases to around 10°, 
alleviating the super-Eddington tension. Many ultraluminous X-ray sourcedI9 have been ob- 
served to have Eddington luminosity ratios above 100, but still, the mechanism of super-Eddington 
accretion remains an open question! If the accretion origin is valid for SN 2022jli, it provides a 


new environment and opportunity to study super-Eddington accretion. 


Under the scenario that accretion powers the excess emission of SN 2022jli, it is intriguing to 
consider the possibility of one pre-existing accreting compact object in the system. Regardless of 
the origin of the compact object, one strong constraint on the accretion is that we need hydrogen- 
rich material to fuel the accretion to explain the hydrogen emission and its evolution in the late-time 
spectra (see below Method Section § 9 and § 10). However, the non-detection of hydrogen lines 
in the photospheric spectra (see below Method Section § 8) indicates hydrogen-poor ejecta, which 
excludes the supernova ejecta as the main fuel for accretion. This means that accretion of the 
supernova ejecta onto a pre-existing compact companion to the progenitor star can not explain the 
spectroscopic properties of SN 2022jli. Therefore, we suggest a newly formed compact object in 
the supernova explosion and a companion star with a hydrogen-rich envelope is the most likely 


binary system. 


8. Spectral analysis 


Photospheric spectra: Two spectra taken on 2022 May 11 (+6.0 days after discovery) and 2022 
May 24 (+19.2 days after discovery) were used to classify SN 2022j1H#213) these spectra are 
available on TNS. We estimated an ejecta velocity of around 8,200 km s~! from the absorption 
minimum of the identified absorption lines. We compared the +19.2d spectrum of SN 2022jli to 
other supernovae (Extended Data Fig. |2) and found that SN 2022jli resembles normal SNe Ib/c 
well, and it is a genuinely good match to the SE SN 2013ge. Helium absorption was found in 
the early-time spectra of SN 2013ge through careful analysis of both the optical and NIR spectral 
sequence, which likely results from a thin layer of helium remaining at the time of core collapse!) 
Actually, a complete stripping of the He layer from the progenitor stars of SNe Ic is not expected 
in many models, and the contribution of the remaining helium layer to the spectra of SNe Ic has 
been long debated H7] Helium absorption might also exist in early-time spectra of SN 2022jli, 
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but we do not have an extensive enough early-time spectral sequence to explore this fully. We point 
out that the absorption signature of Helium was indeed reported in the NIR spectra of SN 2022jli 
presented in refil 


Nebular spectra: Normal SNe Ib/c start to enter the nebular phase some months after the ex- 
plosion, during which the supernova ejecta become optically thin, and emission lines with little 
continuum emission dominate the spectra. The spectra of SN 2022jli are well sampled at late times 
during the long-lasting gradual decline phase and the fast-decline phase. The late-time spectra of 
SN 2022jli show significant differences with those of normal SNe Ib/c, for example, SN 2013ge 
in Supplementary Information Fig.|5| Before the fast-decline phase, the spectra of SN 2022jli are 
good matches to the other two long-lasting SESNe (SN 2012au and iPTF15dtg), showing promi- 
nent permitted Oxygen emission and iron-plateau, and they look more similar to some SLSNe, 
GRB-SNe than normal SNe Ib/c (Extended data Fig. 2). We identify the emission lines in Extended 
data Fig.|4] Besides the aforementioned iron plateau and Oxygen emission, the other elements con- 
tributing to the plethora of emission lines in SN 2023jli spectra can be identified as Ca, Mg, C, and 
Na. The O 1 A7774 line seems to have two different components manifesting as a relatively narrow 
feature on a broader base. The OI A9263 line appears to be an isolated line without significant 
blending from other emissions. We measured the width of the O 1 9263 line by fitting a simple 
Gaussian profile and obtained the velocities as follows: vpwum(O 1\9263, +203.8 d) = 2660 + 90 
km s~!, and vrwym (O 1\9263, +260.9 d) = 2720 +50 kms~'. We also measured the width of the 
narrow component of the O 17774 line, and obtained vpwy(O IA7774, +203.8 d) = 2580 + 200 
kms~!, and vpwum(O 1\7774, +260.9d) = 2480 + 30 kms~!. Before transitioning to the fast- 
decline phase (K +270 days), SN 2022jli shows prominent emission around 6500 A, which differs 
from all the other comparison objects (see Methods § 9 for discussion on this feature). 


Accompanying the fast photometric evolution of SN 2022jli from the gradual decline to the 
rapid drop in luminosity, significant spectral evolution was also noticed (Extended Data Fig. [5p. 
One conspicuous change is the disappearance of those narrow Oxygen lines, including O 1 \6158, 
O1 A7774, O1 48446, and O 1 \9263, and probably also [O 11] A\7320,7330 after the supernova 
enters the fast-declining phase (> 270 days after discovery). The narrow features of the Ca 1 NIR 
triplets also disappeared. In the meantime, the [O I] 4A6300,6363 lines emerge. We measured an 
[O 1] \6300 line width of vrwnm(+280d) = 5520 + 190 kms~'. Theoretically, the fast change 
of OI emission from permitted emission to forbidden emission indicates a dramatic change of 
density or temperature in the ejecta. In the case of SN 2022jli, the transition happened rather 
quickly, during which the ejecta density is not expected to change much. Therefore, the most 
plausible explanation for the spectral change is a decrease in temperature due to the quenching of 
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the central energy source, which is consistent with the contemporary fast decline of the light curve. 
The narrower widths of the disappearing lines corroborate the idea that they were emitted from the 
inner parts of the ejecta with lower velocities. 


Theoretical calculations of nebular phase spectra of SESNe powered by both Nickel decay 
and a central engine are rard O Detailed modeling of SN 2022jli spectra considering both 
energy sources and their temporal evolution might shed more light on the explosion mechanism, 


but is beyond the scope of this paper. 


9. Hydrogen emission in the late-time spectra The spectra of SN 2022jli show unique strong 
emission around 6500 Å (Extended Data Fig. 2). The emission seems to comprise two compo- 
nents, with one narrow feature sitting on top of a broad component. The wavelength range of 
interest in the supernova spectrum is contaminated by host galaxy emission from nearby H HI re- 
gions, but we confirm that the narrow feature does not come from any artifact due to host galaxy 
line contamination (Supplementary Information Fig. (6). We attribute the narrow emission feature 
as Ha emission, and the co-evolution with the emission at the wavelength of H8 supports this 
identification (Supplementary Information Fig.[7). The narrow feature shows back-and-forth shifts 
in wavelength, i.e., moving to shorter and longer wavelengths around the rest-wavelength of Ha 
(bottom left panel of Fig. B). The shifting behavior shows a cyclical pattern, consistent with the 
12.4-day period as derived for the undulation period in the optical light curves. The luminosity of 
the narrow Ha line closely follows the total luminosity of the supernova as Lya = 0.004 x Lpol 
(Extended Data Fig. (10). There is a transition phase in the bolometric luminosity, an overshooting 
before it drops, during which the Ha luminosity also follows the bolometric luminosity closely. 
The narrow feature disappears after the SN enters the fast-decline phase (Extended Data Fig. 5). 
The tight connection between the narrow Ha feature and the light curve of the supernova, both 
the periodicity and the luminosity, implies that the energy that powers the extra emission of the 


supernova is also responsible for the Ha emission. 


The region producing the periodic undulation in the light curve has been limited to a rela- 
tively small size in the center of the ejecta, which implies the emission of Ha also comes from 
the center of the ejecta. The hydrogen material that gives rise to the narrow Ha and H8 emis- 
sion most likely comes from the envelope of the companion star that is accreted onto the newly 
formed compact object. Ha and Hp emission lines have been commonly observed in binary sys- 
tems with accretion disks!2"23] which could serve as an analogy in the low accretion rate regime 
to understand the emission mechanism and structure of the Ha in SN 2022jli. 


After the narrow Ha emission vanishes, there is still a prominent emission around 6500 A (Ex- 
tended Data Fig. [5). Similar emission has been observed in other SESNe!2#26. mainly SNe IIb 
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or Ib. The origin of such emissions has been debated. For example, the emission has been ex- 
plained as either [N 11]A\6548, 6583 or Ha [O] In SNe IIb, the hydrogen could be left- 
overs in the outer layer, and the ejecta-wind interaction could be a possible power source for 
the ionization. However, Dessart et al. (2021)84 argue that the hydrogen envelopes of Type IIb 
SNe are too small and dilute to produce any noticeable Ha emission or absorption after ~150 
days. We performed spectral decomposition of the spectra around 6400 A. The spectra can be well 
decomposed into four Gaussian emission profiles, among which two are emission lines of [O 1] 
A6300, 6363 (Supplementary Information Fig. (8), but the derived velocities for the different lines 
are not consistent for [O I] and the [N 11] or Ha components. We get voy = —890 km s~1(+271d), 
voy = —480kms~!(+275d), voy = —310kms~!(+280d); unm = —3770kms~!(+271d), 
yn = —2690 kms~!(+275d), vnn = —2200kms~1(+280d); vaa = —3260kms~1(+271d), 
VHa = —2170kms~!(+275d), vga = —1680 kms~'(+280d). It is difficult to explain the velocity 
difference if the emission is dominated by [NII]. Another possible origin of the emission could 
be hydrogen stripped from the companion star by the supernova ejecta224| Late-time Ha emis- 
sion with the potential origin of stripped hydrogen from the companion star has been observed 
in several Type Ia supernovae!351137] However, the observed Ha profiles in those SN Ia spectra 
have much narrower line widths around 1000 kms~! while the component on the red side of [O I] 
6300, 6363 in SN 2022jli has a width of FWHM ~ 10,000kms~! which is twice the width 
of the [O 1] 46300 line. The hydrogen origin from the companion star has also been proposed to 
explain the nebular Ha emission in the Type Ic SLSN iPTF13ehe 38139) for which other works 


argue for a hydrogen origin from hydrogen-rich CSM produced in mass loss before the supernova 
explosion 138! 


10. Evolution of the accretion-powered Ha emission In the accretion-powered supernova sce- 
nario, the compact remnant and the companion star are bound in an eccentric orbit, where the 
hydrogen-rich material is accreted from the envelope of the companion star to the newborn com- 
pact remnant every time the compact remnant passes through the pericenter of the orbit. Now we 
consider whether the orbital motion of such an eccentric orbit can explain the observed velocity 
shift. The light-of-sight velocity of the compact remnant can be written as 


V, = V; + K(cos(w + f) + ecos(w)), (2) 


where Vz is the proper motion velocity of the binary system barycenter, K is the characteristic 
velocity amplitude, w is the argument of periapse, f is the true anomaly, and e is the orbital 
eccentricity. We adopted the same nomenclature for the orbital elements as used by ref (40) The 


velocity amplitude can be written as 


mg nasint M2 \9 Me + Mo 2 sini x 
= = 214 ; cng 3 
me + M2 V1— e CeMo ‘I4Mo+15Mo ye (3) 
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Qn 
P 9 
major axis of the elliptical orbit, and 2 is the inclination of the orbit. We can see the velocity 


where me is the compact remnant mass, Mə is the companion star mass, n = a is the semi- 


amplitude is highly degenerate with respect to the binary masses and the orbital eccentricity. 


We extracted the accretion-powered Ha by simply subtracting the emission by a pseudo con- 
tinuum. The continuum model is constructed by fitting a linear model to the arbitrarily selected 
continuum region on both sides of the emission feature. We estimated the Ha velocities by measur- 
ing the flux-weighted centroid of the emission feature without considering the nontrivial velocity 
structure. We estimated the uncertainties for the velocity by considering both the flux uncertainty 
and the systematic uncertainty due to the different choices on the continuum region. The uncer- 
tainty budget did not include the systematic uncertainty that would be introduced by the unknown 
nature of the underlying baseline for the Halpha feature. The result is shown in Extended Data 
Fig. Only the last bump period was well sampled with decent spectral resolutions. We caution 
the unusual rebrightening of the last bump (Method Section § 6), which might introduce other con- 
tributions to Ha velocity besides the orbital motion. Generally speaking, the Keplerian orbit with 
high eccentricity can roughly account for the velocity evolution trend, as shown by four models 
with representative parameters. In this exercise, we simply subtracted one arbitrarily chosen linear 
“continuum” to get the Ha emission without considering the potential contamination from other 
lines. The spectral data analyzed here span more than 120 days, and the spectral evolution could 
cause problems when we combine data from different periods with significant separation in time. 
This might explain the obvious outlier from the earliest spectrum. 


In this heuristic experiment, we assume the orbital velocity is the major component. We 
need to point out that other possible kinematic processes may complicate the velocity evolution. 
These complications include the velocity structure of the accretion disk, the mass flow from the 
companion star to the compact object, and the outflow of winds. The potential eclipse of the 
accretion disk by the companion star and the emission from the companion envelope could also 
complicate the velocity structure of the Ha emission line. We stress that the velocity measurement 
with the current method under simplified assumptions is helpful in understanding the property of 
the binary system qualitatively. We caution the readers to refrain from deriving the exact orbital 


parameters from the current analysis. 


11. X-ray and radio observations: 


Chandra We obtained three epochs of X-ray observation with the Advanced CCD Imaging Spec- 
trometer (ACIS) of the Chandra X-ray Observatory from 2023 January 17 to 2023 January 25 
(+257.4, +261.8, and +265.8 days after discovery) under an approved Director Discretionary Time 
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Proposal (PI: Chen). Each epoch has an exposure time of 10.06 ks. The Chandra ACIS-S data 
were reduced with the CIAO! software package (v4.14) and relevant calibration files (CALDB 
version 4.10.2), applying standard filtering criteria. SN 2022jli was not detected in any of the three 
epochs. We measured the count rate within a 5” radius aperture at the supernova position and ob- 
tained absorbed flux upper limits in 95% confidence interval of 5.27 x 107! erg s7! cm~? (4257.4 
days), 5.32 x 1071 erg s7! cm~? (+261.8 days), and 5.27 x 10715 erg s™t cm? (+265.8 days). We 
merged the three epochs and obtained an absorbed flux upper limit of f,1(0.5 — 7keV;95%) = 


1.77 x10 “ergs “em, 


The Galactic neutral hydrogen column density in the direction of the SN is Ny,,, = 3.49 x 
10% cm~? (ref), From our optical spectra, we estimate E(B — V )pos 


Galactic dust-to-gas ratid“], Ny(cm7?) = 2.21 x 10?4Ay(mag), the extinction value corresponds 


= 0.25 mag. Assuming a 


to an intrinsic neutral hydrogen column density of Nu, ~ 1.7 x 10?! cm~?. For an assumed 
simple power-law spectral model with spectral photon index l = 2, we find an unabsorbed flux 
limit in 95% confidence interval of 2.18 x 10~'° erg s7! cm~? (0.5-7 keV). At the distance of 22.5 


Mpc, this flux translates into a luminosity of 1.32 x 10°° erg s71. 


The supernova ejecta may be optically thick to soft X-rays for decades. The opacities of X- 
ray flux absorption are dominantly due to photoelectric absorption (photon-ionization and photon- 
excitation) below around 100 keV. The X-ray optical depth at a given epoch can be described as 


T(E) ~ me where Mej is the supernova ejecta mass, vej is the ejecta velocity, and «aş is the 
bound-free opacity which is largely determined by the ejecta abundance. Alp et al calculated 
a typical optical depth of T(t, E) ~ 100¢,;?E~? for core-collapse supernovae, where t4 is the 
time since the explosion in units of 10000 days, and E is the energy in units of keV. In the case 
of SN 2022jli, we have 7(+250d,1keV) = 1.6 x 10°, 7(+250d,10keV) = 1.6 x 10°. The 
supernova progenitor in their calculation has an ejecta mass of around 10 Mo. Even if SN 2022jli 


has a smaller ejecta mass, the X-ray optical depth in the Chandra energy band is still very high. 


In the context of interacting supernovae, our upper limits constrain the X-ray luminosity to 
be more than three orders of magnitude lower than the optical/NIR luminosity at the same epoch. 
Such a ratio of X-ray to optical/NIR luminosity puts a strong constraint on the interaction contribu- 
tion to the bolometric luminosity. Higher X-ray to optical/NIR luminosity have been observed in 
supernovae with strong interaction, for example, Type IIn SN 2010j122145) Type Ibn SN 2006j¢!44) 
and Type IIn SN 2006jd 47148) where strong X-ray emission from ejecta-CSM interaction is de- 
tected by the Chandra observatory. 
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NuSTAR SN 2022jli was also observed with the Nuclear Spectroscopic Telescope Array (NuS- 
TAR 4), NuSTAR has two coaligned X-ray telescopes, with corresponding focal plane modules 
FPMA and FPMB. The observations were conducted in three epochs spanning from +227 to +237 
days after discovery. The first epoch (ID 90801535002) was on 2022 December 18, with an ef- 
fective exposure time of 20,353 seconds (FPMA) and 20,142 seconds (FPMB). The second epoch 
(ID 90801535004) was on 2022 December 23, with an effective exposure time of 20,231 seconds 
(FPMA) and 20,008 seconds (FPMB). The third epoch (ID 90801535006) was on 2022 Decem- 
ber 27, with an exposure time of 20,257 seconds (FPMA) and 20,036 seconds (FPMB). The data 
were reduced using HEASoft v.6.31 and the NuSTAR Data Analysis Software (NuSTARDAS) 
v.2.1.2, in particular the nupipeline (version 0.4.9) and nuproducts routines. No 
source was detected at the supernova position. We focused our analysis on the 30-60 keV energy 
range because the probability of detecting the supernova is higher in the hard X-ray energy range 
with lower optical depth. We estimated an upper limit of 4.5 x 1074 countss~! by calculating 
3 x VBriot/ texp Where Bot is the total count in a circular aperture with 50” radius at the position of 
SN 2022jli from both epochs and both instruments and texp = 121 ks is the total effective exposure 
time. Assuming a power-law model with photon index I’ = 2, the above count rate upper limit cor- 
responds to an unabsorbed flux of f,1(30 — 60 keV) = 4.1 x 107! ergs”! cm~? or an upper limit 
of 2.5 x 10% erg s~* at a luminosity distance of 22.5 Mpc. 


ATCA We observed SN 2022jli with the Australian Telescope Compact Array (ATCA in the 
C/X-band (3.9-11.0 GHz) on MJD 59916 for 2 hours and on MJD 59931 for 5 hours using the 
extended 6C array configuration (CX517). For both observations, we used the ATCA calibrator 
source PKS 1934-638 as the primary flux calibrator and the calibrator source PKS 0003-066 for 
phase calibration scans. We used the Mi riad software!="|to reduce and image the data from these 
observations. We used mfclean with Briggs weighting and a robust parameter of 0.0, with 


multi-frequency synthesis (mf s) deconvolution. 


Neither observation resulted in a clear detection of SN 2022jli in the observed frequency 
range. We use imfit in C-band (5.5GHz) to extract 50 non-detection limits of < 0.176 mJy 
beam™! and < 0.074mJy beam™!, respectively on MJD 59916 and MJD 59931. At the X-band 
9.0 GHz central observing frequency we obtain a 5o limits of < 0.057 mJy beam~! and < 0.055 mJy 
beam™!, respectively on MJD 59916 and MJD 59931. 


12. Light travel time and diffusion time We can use the time scale of the 12.4-day period and 
the time scale of around three days in the rising part of the undulation profile to constrain the size 
and location of the varying energy source. If the process that generates the energy powering the 
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periodically undulating light curve happens outside of the supernova ejecta, for example, the ejecta 
and CSM interaction (ECI), the finite size of the ejecta can cause different light travel time from 
different parts of the ejecta, i.e., the emission from the closer side to the observer arrive earlier. In 
the case of ECI, if we consider a simple model where the optically thin ejecta expand homologously 
and the CSM has a spherical shell-like structure, the light travel time difference Hia the earliest 


and latest arrival of emission from the same sphere is At ~ z = 10 x =—days. Such 


104 mae 1 mi 
a phase-dependent light travel time difference will smear out any potential periodic signals. 

If the energy source is located inside the ejecta, any sudden change in the input luminosity is 
subject to delay and smear effects by the ejecta that the emission goes through before reaching the 
observer. The short rising time of around 3 days of the bump requires the diffusion time to be less 
than that. The diffusion time through spherical ejecta with opacity « from a radius r to the outer 
edge R is 


fk 
tag(r, R) = “| p(r)rdr. (4) 


If adopting ejecta density profile p(r) = = =a (ref!52), we get 


M 
ta(r, R) = a (1-5) . (5) 


For ejecta in homologous expansion, we have R = vt where v is the velocity of the outer 
edge of the ejecta, and t is the time since the supernova explosion. We get the diffusion time from 
the center of the ejecta 


E 3's g ( ° yo t Vy (6) 
as a Mo 0.07 cm? g7! \104 kms7! 50 day A 


If we assume the discovery time of SN 2022jli is around the first peak and assume a typical rising 
time of 15 days for SNe Ic, the second peak is around 70 days after the explosion. We note that 
the discovery time could be later than the first peak because we missed the rising part of the first 
peak. Therefore, at the second peak, adopting ejecta velocity of 8,200 kms, the requirement of 
diffusion time less than 3 days is roughly consistent with the energy source residing in the center 
of the ejecta if the ejecta mass Mej S 1.7 Mo. 


13. Supernova explosion in a binary system Two mainstream scenarios, the stellar-wind origin, 
and the binary interaction origin, have been proposed to explain envelope stripping for the pro- 
genitor of SESNe. Observational evidence has been found in both directiond'=!54 A supernova 
explosion in a compact binary system, i.e., a binary system with a small separation between two 


stars, significantly impacts the companion star and the further evolution of the binary systen 5 
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The existence of the companion star, in return, might have profound effects on the manifestation 


of the observed supernova. 


Assuming the mass loss in the supernova explosion is instantaneous, the equations relat- 
ing the pre-explosion and post-explosion orbital parameters (in an instantaneous reference frame 


centered on the companion star right at the time of explosion) are 155) 


Vi GMM, _ GMM, cn 
2 Qi 2aş 


and 


ugay fw? + (Vi + w)? = upy G(M. + Ma)ag(1 — €) (8) 


where M, is the companion mass, Me is compact remnant mass, pp = MeMa/(Me + Mə) is 
the reduced mass of the system after explosion, V; is the relative velocity before explosion, V+ is 
the relative velocity after explosion, a; is the semi-major radius of the orbit before explosion, ap is 
the semi-major radius of the orbit after explosion, Ù = (w,, Wy, wz) is the kick velocity. 


The above equations result in the following evolution of orbital parameters: 


2 2 2 -1 
ar w, ee + (Vi + wy) 
Jf _ |9_ 
ay | x( V? a 
and 5 (v, j2 
Lee" = eat 22 : d 10 
Xag ( V? ) ca 


where x = (Mı + M2)/(M. + M2) > 1. M; is the mass of the primary star that exploded. When 
the explosion is spherically symmetric, the compact remnant gains no kick velocity, i.e., Vp = V;. 
If the binary system remains bound, we get the simplified results of the orbital parameters of 


af 1 


baad ag = 11 
a Doe (11) 
and M.—-M 
1 Cc 
eee 12 
e M. 4 M, (12) 


The above results tell that the supernova explosion will widen the binary orbit and introduce ec- 
centricity to the system even when the explosion is symmetric. 


In the case of SN 2022jli, the post-explosion binary system has 


(13) 
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where P = 12.4 days. Before the supernova explosion, the separation between the primary and 
secondary stars is smaller, which can easily achieve Be < 10, where R» is the radius of the com- 
panion star, for a large parameter space of Mı, Mə, and Me, in which regime the ejecta impact has 
a significant influence on the companion! 3311561158] The momentum transfer and energy injection 
will result in an impact velocity to the companion, strip material from the companion, and bloat 
the envelope of the companion. The amount of removed stellar mass, the resulting impact velocity, 
and the companion’s reaction to the impact in radius/temperature/luminosity strongly depend on 
the binary separation and explosion energy. For a main-sequence companion star, the radius can 
easily increase by one or two orders of magnitude 331158] for a typical explosion energy of 10°! erg 
and 7 < 10. The fluffy material of the bloated envelope makes it easier to get accreted onto 
the supernova remnant, which could form an accretion disk and provide extra energy to power the 
supernova”, The accretion rate is not expected to be constant since the companion radius evolves. 
Once the deposited energy fully radiates away, the companion starts to contract by releasing their 
gravitational energy and eventually resumes their original state before the SN explosion. 


Due to accretion from the companion star to the compact remnant with likely wind and jet 
production in the meantime, the evolution of the surviving binary system is subject to the mass 
loss and mass transfer process. The detailed orbital evolution of the binary system is beyond the 
scope of this work. But we can argue qualitatively if any significant orbital evolution is expected 
during the ~ 300-day evolution of SN 2022jli. If we assume non-conservative, i.e., M, Jab # 0 
where M is the total mass change of the system and Joab is the orbital angular momentum, delta- 
function mass loss/transfer through Roche Lobe Overflow (RLOF), the secular relative change of 
the orbital semimajor axis has (2) ~ ea with an order-of-magnitude estimate 52, where M, is 
the mass of the accretor and M; is the mass of the donor. With dM) ~ Mace = 107° A Mo 
for SN 2022jli (Method Section § 7), and assuming Mı = 5Mo, we can tell the relative change of 
semi-major axis during the first ~ 300 days of SN 2022jli is very small which is consistent with 
what we found in Method Section § 5. 


A surviving companion star has been important evidence for the binary origin of some 
SESNe. Such companion stars might have been observed in SN 1993} (SN IIb), SN 2001; 
(SN IIb), SN 2006jd!8!1@2l (SN Ibn), and SN 201 1dh!8I!@4l (SN IIb). SN 2022jli provides a good 
candidate to search for the surviving companion star of a Type Ib/c supernova. 


Detailed numerical simulations have shown that a supernova explosion can often unbind the 
secondary star in the binary system 651/66] In fact, re?l found that around 90% of the pre-explosion 
binary systems become unbound after the core-collapse explosion, and the rarity of a bound binary 
system surviving a supernova explosion has been observationally confirmed by dedicated searches 
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for the companion star of the compact remnant in young supernova remnants!©71'68| and search for 
companion stars of magnetars=4) This may be the reason why we did not observe similar properties 
of SN 2022jli in other SESNe. 


14. Possible explanations for the late-time rapid decline One interesting phenomenon of SN 2022jli 
is the rapid decline around 260 days after the discovery. One possible explanation is that the rapid 
decline happens when the accretion rate drops significantly after the bloated envelope of the com- 
panion star shrinks back (Method Section § 13). The other possible explanations involve a dramatic 
change in the binary orbit, for example, the compact object plunges into the companion star B 
which prevents further mass accretion. For the latter scenario to happen, the compact object is ex- 
pected to penetrate the envelope of the companion star and lose significant angular momentum and 
mechanical energy to shrink the distance to the companion star at periapse. For an eccentric orbit, 
the periapsis distance is rperi = a(1 — e). The period of the binary system in SN 2022jli, there- 
fore the semimajor axis, did not show any significant evolution before the rapid decline (Method 
Section § 5). The mechanical energy of a binary system is E = -A and the absolute value 
of the orbital momentum is Jop = p,/GMa(1 — e?) where M = Mı + M; and u = mMm, If 
we simply assume an external force working on the compact object to slow the velocity instan- 
taneously at the periapse, the mechanical energy of the binary system is more susceptible to the 
external force than the angular momentum. It is unlikely that the orbital eccentricity changes sig- 
nificantly while the semimajor axis remains constant. In the first explanation, the binary orbital 
parameter undergoes secular evolution, which barely changes in the comparatively short duration 
of the supernova. The accretion rate suffers a fast drop, but the periodicity of the orbital undulation 
shall persist if the compact object continues to accrete mass from the companion. We inspect the 
late-time evolution of the pseudo-bolometric light curve of SN 2022jli to see if the 12.4-day peri- 
odicity continues. During the fast decline phase, there is a hint of a small bump in the light curve 
as indicated by the vertical dashed line in panel b of Extended Data Fig. [6] which is consistent 
with 12.4-day separation to the peak of the last clear bump before the rapid decline. This might 
be a hint of evidence for the continuous accretion. The prominent accretion rate variation before 
the rapid-decline phase might change the disk properties significantly. We suspect this could trig- 
ger some disk instability”, which might be responsible for the enhanced luminosity before the 
rapid-decline phase, as reported in Method Section § 6. 


15. Fermi-LAT detection of y-ray emission from the direction of SN 2022jli The Large Area 
Telescope (LAT) on the Fermi Gamma-Ray Space Telescop! has been surveying the entire sky 
since 2008. LAT has a large field of view of ~ 60°, enabling it to scan the sky in about 3 hours. 
We queried LAT data within a 10° radius of SN 2022jli with photon energies between 100 MeV 
to 300 GeV observed in the past 14.5 years between 2008 September 1 and 2023 March 1 us- 
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ing the Fermi Science Support Center data servel We filtered the photons using source-type 


events (evclass=128) with the most stringent cuts on the data quality (DATA_QUAL==1 && 
LAT_CONF1IG==1), reconstructed both in the front and the back of the detector (evt ype=3) and 


with a maximum zenith angle of 90°. 


New y-ray source detection We used the standard binned likelihood analysis method gt Like 
to analyze the data. The sources in the incremental Fermi Large Area Telescope Fourth Source Cat- 
alog (4FGL-DR3; 27), together with the diffuse Galactic and isotropic backgrounds!lg1 l_iem_v07.fits 
and iso_P8R3_SOURCE_V3_v1.txt, are included in the model. We tried different binning 
strategies in time and energy to search for new y-ray sources. A significant new source was de- 
tected in data observed after the supernova explosion from 2022 May 1 to 2023 March 1, as shown 
in the test statistic (TS) map in the top left panel of Extended Data Fig.|7| The energy light curve 
of the new source in the 100 MeV — 300 GeV energy band is shown in the bottom left panel of Ex- 
tended Data Fig.|7| The source was detected in two bins corresponding to September and October 
2022 (TS=20.0) and November and December 2022 (TS=29.6). We got the spectral energy distri- 
bution of the new y-ray source in three different time windows: September and October 2022 (top 
left panel of Supplementary Information Fig. (9); November and December 2022 (top right panel of 
Supplementary Information Fig. (9), and the above two combined (bottom panel of Supplementary 
Information Fig. (9). The new source is most significantly detected in the 1—3 GeV energy band. 
We also extracted the 1—3 GeV light curve at the SN position, and the result is shown in the bottom 
right panel of Extended Data Fig. 


We performed likelihood modeling with the 1-3 GeV data observed from 2022 November 
1 to 2023 January 1. The TS map is shown in the top right panel of Extended Data Fig.|7} and the 
localization of the new source from the above modeling is shown in Extended Data Fig. |8} The 
best localization of the new y-ray source is RA= 8.620(+0.084)°, Dec= —8.425(+0.077)° with 
uncertainty given by the 68% confidence interval. The Fermi-LAT source 4FGL JO035.8-0837 at 
RA=8.958(+0.060)°, Dec= —8.632(+0.048)° is 0.33° East, 0.21° South to the newly detected 
source. They are two different sources spatially separated from each other, but we note that they 


could contaminate each other when measuring their fluxes. 


Poisson noise matched filter method We also used the Poisson noise matched filter to corrob- 
orate the detection of the GeV photons associated with SN 2022jli. Since we are dealing with the 


problem of detecting sources embedded in low-number-count Poisson noise, the optimal matched 
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filter presented in ref" was used. We built the Poisson-noise optimal filter as 
p p 


F 
Pra =In (1+ 5P) ; (14) 


where B is the expectancy for the background level, F is the unknown flux of the source we 
would like to detect, and P is the PSF of LAT (the observed photon distribution of a point source 
with LAT). We divided the 14.5-year data into two-month bins, and we then divided the pho- 
tons into seven energy bands (0.1 —0.3 GeV, 0.3 — 1.0 GeV, 1 —3 GeV, 3 — 10 GeV, 10 -30 GeV, 30 — 
100 GeV, 100-300 GeV) in each temporal bin. We generated count maps for each temporal and 
energy bin. Following the method in ref!) the log-likelihood difference image used to detect the 
source is 


S=M@ Pr (15) 


where M denotes the measured data, i.e., the count map, ® denotes convolution, + denotes coor- 
dinates reversal. The PSF of Fermi-LAT observation is a function of an incident photon’s energy, 
the inclination angle, and the event class. In practice, we built the PSF by generating energy- 
dependent PSF data with gt psf for the two-month observation and then simply averaged the PSF 
of different energies. 


The result shows a significant source detection at the SN position in the 1—3 GeV band in 
November and December 2022, as shown in Supplementary Information Fig. which is consis- 
tent with the result obtained with gt like. The other Fermi-LAT source near SN 2022jli, 4FGL 
J0035.8-0837 has also been clearly revealed in the 1—3 GeV band in March and April 2020 (Sup- 
plementary Information Fig.|11), and also in the 3—10 GeV band in March and April 2020, and 
July and August 2020 (Supplementary Information Fig. [15). Emission of 4FGL J0035.8-0837 
might also contribute to the photons seen in July and August 2013 (Supplementary Information 


Fig. [14). 


Evidence for periodicity of the new y-ray source We studied the temporal distribution of the 
photons associated with the new y-ray source. There are eleven 1—3 GeV photons within the half 
containment radius of the averaged PSF from 2022 September | to 2023 January 1. The estimated 
number of background photons within the above temporal and spatial range is two. As shown in 
Extended Data Fig.|8} there is a large avoidance region (50%) in the phase space after folding the 
y-ray photon light curve with a 12.4-day period, as shown by the shaded grey area. This hints 
at a correlation between the y-ray photon arrival time and the optical light-curve bump phase. 
We performed a simple simulation by drawing N photons randomly distributed in four months, 
assuming there is no preference for when the photons arrive with respect to the bump phase. 98.9% 
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(98.0%, 96.5%) of the experiments result in the maximum separation of less than 0.5 for any two 
photons in the phase space for N = 11 (N=10, N=9). We checked and found that the correlation 
of y-ray photons with the 12.4-day period could not be caused by the telescope survey profile, as 
shown in Supplementary Information Fig. 


Association between SN 2022jli and the new y-ray source The new Fermi-LAT y-ray source is 
positionally and temporally consistent with SN 2022jli, which suggests that the new y-ray source 
is associated with the supernova explosion. The correlation between the y-ray photon detection 
time and the 12.4-day periodic undulation of the optical light curve provides further evidence that 
the new y-ray source and SN 2022jli share the same origin. 


Due to the potential periodicity detected in the new y-ray source, the probability of it com- 
ing from some other contamination source, such as a blazar, within the localization area is low. 
Nevertheless, we inspected any potential sources within the 99% localization confidence area 
which might emit 7-rays. First, the nucleus of NGC 157 is a dormant supermassive black hold” 
which does not belong to either flat spectrum radio quasars (FSRQs) or BL Lac objects (BL 
Lacs). The only blazar candidate within the 99% localization confidence area is NVSS J003456- 
082820, as shown in Extended Data Fig. NVSS J003456-082820 was listed as a blazar can- 
didate in the Blazar Radio and Optical Survey (BROS catalog based on its spectral index 
between the radio bands of 0.15 GHz and 1.4 GHz and its compactness. We did not detect 
any optical flare from NVSS J003456-082820 in ZTF data but noted that around 20% of the 
orphan y-ray flares happen without an accompanying optical flar! The true nature of NVSS 
J003456-082820 requires further confirmation from multifrequency observations. We studied the 
expectancy of one BROS catalog object falling in the localization area of the new y-ray source 
by calculating the chance coincidence probability of observing NVSS J003456-082820. We cal- 
culated Py, = 1 — exp[—7(RZ + 40?) pa] where Ro is the angular distance between the new y-ray 
source location and NVSS J003456-082820, and ø, is the 68% localization uncertainties radius of 
the new y-ray source, and p4 is the surface density of BROS sources. The result of Pen œ 0.5 tells 
that having one BROS object in the observed field of the new y-ray source is not surprising. 


Origin of the SNe-associated y-ray emission Several physical processes could generate cos- 
mic y rays. These processes include radioactive decay, particle-particle collisions, acceleration 
of charged particles (bremsstrahlung or synchrotron radiation), inverse Compton scattering, and 
matter-antimatter annihilation. One notable category of y-ray sources is Gamma-ray Bursts (GRBs), 
lasting from ten milliseconds to several hours, which is different from the gamma-ray emission we 


observed in SN 2022jli with a month time scale. Gamma-ray sources with radioactive decay origin 
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are usually associated with nucleosynthesis. For example, y-ray emissions from °°Co decay at 847 
and 1248 keV have been detected in some very nearby supernovae such as SN 1987A and 
SN 2014J "80 The GeV emission in SN 2022jli apparently has different origins. 


Significant efforts have been devoted to searching for other types of y-ray emission asso- 
ciated with supernova explosion 1/8184] These efforts are motivated by different mechanisms 
that could happen in the supernova environment producing y-ray emissions, for example, proton- 
proton collisions in supernova ejecta and CSM interaction == and inverse Compton scattering in 
the young neutron-star wind nebula!l/82| However, the detected y-ray sources associated with 
supernovae are limited. The only tentative detections of such objects are reported for the Type II 
supernova iPTF14hls == and Type IHP SN 2004dj 185] This is the first time significant GeV y-ray 
emission has been detected from a SESN. We have shown that ECI or a magnetar cannot explain 
the observed optical light curves of SN 2022jli, which disfavors y-ray generation from such pro- 
cesses. The accretion of hydrogen-rich mass from the companion star is expected to happen in 
the center of the supernova ejecta. Potential y-ray emission in accretion-powered supernovae has 
barely been discussed in the literature, but y-ray emission generated in accretion-related processes 
in microquasars has been commonly observed, for example, in Cygnus X-3 [86[187] The ultrahigh 
accretion rate required to power SN 2022jli is susceptible to launching extreme outflow winds or 
an accretion jet. We speculate that the observed y-ray photons in SN 2022jli could be produced by 
such a wind or a jet analogous to the observed y-ray emission in some X-ray binaries“ 
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Data availability 


Photometry and spectroscopy of SN 2022jli will be made available via the WISeREP public database. 
Facilities that make all their data available in public archives promptly or after a proprietary period 
include Palomar 48-inch/ZTF, VLT/X-Shooter, NuSTAR, Chandra X-ray observatory, and Fermi 
Gamma-ray Space Telescope. Data from ATLAS, Gaia, ASAS-SN, and KKO were obtained from 


public sources. 


Code availability 


The code and data used to perform the analysis and produce the figures for this paper are available 


in a public GitHub repository (https: //github.com/AtomyChan/SN2022j1i). 
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Extended Data Figure 1: SN 2022jli and the host galaxy NGC 157. The background image 
shows the 3-color (Y, H, K bands) image of NGC 157 taken before the supernova explosion with 
the HAWK-I instrument on ESO’s Very Large Telescope (VLT) at the Paranal Observatory in Chile 
(Credit: ESO). The red plus symbol in the figure indicates the position of SN 2022jli. The inset 
panels: (a) zoom-in view around SN 2022jli showing the nearby environment in NIR; (b) zoom-in 
view around SN 2022jli on an r-band image taken with Magellan/IMACS on 2022 December 15; 
(c) NaI absorption lines from the host galaxy; (d) the narrow diffuse interstellar band absorption 
(DIB6283) from the host galaxy. Panel (a) and (b) share the same field of view size. In panels (c) 
and (d), the red spectrum shows the averaged spectrum of three X-Shooter spectra and one IMACS 
spectrum, whereas the individual spectra are shown in grey in the background. 
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Extended Data Figure 2: Spectral evolution and comparison of SN 2022jli. (a) Photospheric 
spectra of SN 2022jli. The +6d spectrum and the +19d spectrum are shown in black and red, re- 
spectively. The identified absorption lines of Fe II, OI, and Ca II are indicated with the arrows. 
A blue-shifted velocity of 8, 200km s™t has been applied to all the lines. (b) The +19d spec- 
trum of SN 2022jli compared with photospheric spectra of other SNe. The comparison objects 
include: normal SNe Ic (SN 2004g K3] SN 2007g 9, SN Ib/c (SN 2013ge!4h, “Hypernovae”’ 
(SN 1997dq PEED, Hydrogen-poor SLSNe (SN 2015br2!4), and other long-lasting peculiar 
SESNe (SN 2012au", iPTF15dtg "5). (c) Late-time spectra of SN 2022jli and other compari- 
son SNe. The comparison objects are the same as in panel (b) except that we have added another 
Hydrogen-poor SLSN SN 2007bi © The verticag slashed lines mark the wavelength of Ha and H2. 
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Extended Data Figure 3: Periodicity analysis of light curves in the individual bands. Left: the 
Lomb-Scargle power spectrum of the light curve (blue) and the observation window function (yel- 
low). The zoom-in panel shows the region of interest around the true periodic signal. The hori- 
zontal lines indicate the false alarm probability (FAP) levels of 1073, 1076, and 107° from top to 
bottom. Right: the phase-folded subtracted and normalized light curve. The period adopted for the 
folding corresponds to the peak power in the zoom-in panel on the left. 
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Extended Data Figure 4: Optical and NIR spectra of SN 2022jli around +210 days after discovery 
(optical spectrum taken with NOT/ALFOSC on 2022 November 25 and NIR spectrum taken with 
Magellan/FIRE on 2022 December 15). The NIR spectrum has been scaled to match the flux of the 
optical spectrum in the overlapping region. The prominent features and the corresponding ions (or 
molecules) that likely contributed to the emission lines are marked. The spectrum in blue shows 
the mock spectrum of Fe II, Fe IN and Ni II emission, and the vertical lines indicate emission from 
individual transitions which give relatively strong emission. The absolute strength of the mock 
spectrum is arbitrary, and the relative strength between different transitions for the same ion was 
calculated for a temperature of 104 K, and an electron density of 10’cm~°. 
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Extended Data Figure 5: Spectral evolution of SN 2022jli before and during the fast decline 
phase of the light curve. The top left panel shows the bolometric light curve, and the top right 
panel shows the zoom-in view of the fast decline phase of the light curve. The second panel shows 
the six epochs of spectra as indicated in the top right panel. The bottom three panels show the 
zoom-in view of the spectra within different wavelength ranges. All spectra in the bottom three 
panels have been scaled to have the same integrated luminosity between 3800 to 9000 Å as the 
BINOSPEC spectrum taken on 2022 January 12. 
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Extended Data Figure 6: Pseudo-bolometric light curve of SN 2022jli. (a) The points show the 
pseudo-bolometric light curve from 3750 A to 25000 A. These data points share the same colors 
with the other panels to indicate different phases. The black shows the gradual decline phase 
during which the constant relative undulation is detected. The last bump before the fast-declining 
phase is shown in blue, and the fast-declining phase is shown in cyan. The magenta line shows the 
radioactive decay model with 0.15 Mo *°Ni. (b) Zoom-in view of the accretion-powered pseudo- 
bolometric light curve, Laccretion = Luvorr — Lss yi, before and during the fast-declining phase. 
The black solid line shows the linear fit to data between 200 and 260 days after JO=2,459,700. 
The red solid line shows the best-fit undulation model to the data, while the red dashed line shows 
the extrapolation of the undulation model if the SN follows the previous undulations. The vertical 
green lines mark the 12.4-day periods with the left three lines at the minima of the undulation 
profiles. (c) The undulation profiles of the gradual-declining phase and the last bump. (d) The 
relative undulation of the accretion-powered pseudo-bolometric light curve. The red line shows 
the empirical undulation model adopting the same empirical undulation profile as in panel (c). 
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Extended Data Figure 7: Detection and light curve of the new y-ray source. The top panels show 
the Test Statistic (TS) map of the region of interest in the direction of SN 2022jli. (a) The result of 
data observed from 2022 May 1 to 2023 March 1 in the broad energy band (100 MeV — 300 GeV); 
(b) The result of data observed from 2022 November 1 to 2023 January 1 in the narrow energy 
band (1-3 GeV). The sources from the 4FGL-DR3 catalog shown with blue plus symbols have 
been modeled and subtracted from the map. The bottom panels show the y-ray light curves of the 
detected source at the position of SN 2022jli with a bin size of 2 months. (c) The energy range 
of 100 MeV to 300 GeV; (d) The energy range of 1 GeV to 3 GeV. The black points show the 
measured energy flux from the likelihood modeling with the Fermi-LAT analysis tool. The blue 
lines give the upper limit of energy flux within a 95% confidence interval. The red histograms 
show the Test Statistics values on the right axis. The horizontal dashed red line marks TS=9. The 
vertical yellow lines mark the discovery time of SN 2022jli. 
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Extended Data Figure 8: Localization and potential periodicity of the new y-ray source. (a) The 
background mosaic gray pixels show the count map of 1-3 GeV photons detected between 2022 
November 1 and 2023 January 1. The pixel size is 0.125° x 0.125°. The yellow plus symbol 
shows the best-localized position of the new y-ray source, and the surrounding yellow contours 
show the corresponding 68%, 95%, and 99% confidence area. SN 2022jli, the red plus symbol, 
is within the 68% uncertainty region of the detected y-ray source. The blazar candidate NVSS 
J003456-082820 is shown with the blue plus symbol and is within the 95% uncertainty region of 
the new y-ray source. In the central 1.625° x 1.625° field, there is one detected y-ray source from 
the LAT 12-year Source Catalog (4FGL-DR3), 4FGL J0035.8-0837. The black circle has a radius 
of 0.4° corresponding roughly to the 50% containment radius of the averaged PSF over the energy 
range between 1 GeV and 3 GeV. (b): The distribution of 1—3 GeV photons of the new y-ray 
source. The top panel shows the photon energy and detection time. The bottom panel shows the 
distribution after folding the light curve with a period of 12.4 days. The reference time (phase=0) 
corresponds to the minimum of the optical undulation profile. Most of the y-ray photons come 
from the rising phase of the optical bump. The 11 photons are within the half-containment radius 
shown in panel (a). (c): The cumulative distribution of the maximum separation between any two 
photons that would be achieved by drawing N photons randomly distributed in a time range of 120 
days. 61 
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Extended Data Figure 9: Light curve of SN 2022jli compared with those of other supernovae. 
(a) Comparison with Type Ic supernovae dominantly powered by radioactive decay showing clear 
exponential decay tails. All the comparison supernovae have been shifted to have a peak at 14th 
magnitude. The inset panel shows a zoom-in view around the peak light. (b) Comparison with the 
long-lasting SN 2012au (SN Ib), and other supernovae with double-peaked light curves. 
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Extended Data Figure 10: Evolution of the accretion-powered Ha emission. (a) The line lumi- 
nosity of the Ha emission compared with the pseudo-bolometric luminosity. (b) The velocity of 
the Ha emission. The data points share the same color as in panel (a), indicating the phase of the 
corresponding spectrum. The lines show the orbital velocity model with specific orbital parameters 
in the legend. The orbital parameters include compact remnant mass me, companion star mass M2, 
and orbital eccentricity e. The errorbars are 1o confidence intervals. 
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Supplementary Information Figure 1: Light curve of SN 2022jli. The vertical lines on top 
show the epochs of the spectra. The epochs of the NIR spectra are shown in red color, and the 
optical spectra are shown in black. The first two dashed lines indicate the two classification spectra 
obtained from the TNS. The offsets applied to different bands are given in the legend. 
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Supplementary Information Figure 2: Spectra of SN 2022jli taken with P60/SEDM. The phases 
are reported relative to the discovery time of JD = 2459704.67. 
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Supplementary Information Figure 3: Spectral evolution of SN 2022jli. The phases are reported 
relative to the discovery time of JD = 2459704.67. The locations of telluric features are marked 
with a vertical green dashed line. 
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Supplementary Information Figure 4: NIR spectra of SN 2022jli at phase +113, +224, and +258 
days after discovery time. To separate the FIRE spectrum on 2022 December 14 and the X-shooter 
spectrum on 2023 January 18, the X-shooter spectrum has been multiplied by a scale factor of 
0.2. The upper panel shows the transmission curve used to correct the X-Shooter 2023 January 
18 spectrum, which gives an idea of how the telluric absorption affected the observed spectrum in 
different wavelengths, especially the total absorption around 1.4 um (between J and H band) and 
1.85 um (between H and K band) result in large residuals marked by the shaded grey area. 


@) o ae, (c) è — SN 2022jli g 
0] v A 
g 150 ^ SARAS manba aat uy aie: 
= (a) d - 
T 17.54 f PHS, ey $— SN 2013ge V 
D b 3 
2 (b) — 
= 20.04 © "a 
3—$__5_| 
T T T T T T T (d) T 
0 50 100 150 200 250 300 350 400 
Phase 
10} (a) — SN 2022jli +23.9d 


— SN2013ge +9.3d 


a 


~~ 
F 
o 0 
q (b) —— SN 2022jli +67.9d 
E£ 20 4 — SN2013ge +75.3d 
(S) SN2013ge +75.3d x9.3 
T 
7) 
104 
> 
(o) 
2 a ee, Ne 
T 0 
2 (c) —— SN 2022jli +163.8d 
= 7.54 — SN2013ge +160.6d 
D Ha SN2013ge +160.6d x8.6 
c 5.07 
5 
x 25 
= 
WL 


(d) —— SN 2022jli +284.3d 
—— SN2013ge +337.3d 
SN2013ge +337.3d x8.2 


4000 5000 6000 7000 8000 9000 
Rest wavelength [A] 


Supplementary Information Figure 5: Comparison between SN 2022jli and the normal Type 
Ib/c SN 2013ge. Top: g- and r-band light curves of SN 20222jli and V- and r-band light curves of 
SN 2013ge. All the light curves have been corrected for extinction from both the Milky Way and 
the host galaxy, with total extinction of E(B — V) 2022); = 0.289 mag and E(B — V )2013ge = 0.067 
mag. Bottom: (a), (b), (c), and (d) subpanels show the spectra of SN 2022jli, and SN 2013 at 
different phases, as indicated with the arrows in the top panel. In (b), (c), and (d) subpanels, the 
spectra of SN 2013ge shown as dashed lines have been scaled to match the r-band magnitude of 
SN 2022jli at the comparison epoch. 
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Supplementary Information Figure 6: X-Shooter spectra of SN 2022jli in the 5500 — 7500 Å 
region. The bottom panel shows the three X-Shooter spectra after removing artifacts due to over- 
subtraction of host galaxy emission lines and smoothing. The central panel shows the zoom-in view 
of the spectra near the Ha region shown in the box. The top three panels show the wavelength and 
flux calibrated 2D spectra where the over-subtraction of host emission lines ([N 11] AA6548, 6583, 
Ha) is clear. 
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Supplementary Information Figure 7: X-Shooter spectra of SN 2022jli around Ha (in coral) 


and H£ region (in sky blue) showing similar profile and coevolution between these two emission 


features. Left: Normalized spectra in velocity space. The velocity was calculated to the rest-frame 
wavelength of Ha (6562.79 A) and H8 (4861.35 A). The spectra were normalized such that the 
lowest flux is 0 and the highest flux is 1 within the velocity range from — 10,000 km s~! to 10,000 
km s~!. Right: The pseudo-continuum (indicated with the dashed line in the left panels) subtracted 


spectra of the emission feature. 
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Supplementary Information Figure 8: Spectral decomposition of SN 2022jli around 6400 Å 
which are attributed to emission lines of [O 1] \A6300, 6363, Ha, and/or [N 11] AA6548, 6583. The 
rest-frame wavelengths of the relevant lines are indicated by the vertical dashed line in cyan ([O I] 
A6300, 6363), magenta (Ha), and green ([N 1n] \A6548, 6583). The spectra have been normalized 
by the maximum flux within the region shown in the plot and have been shifted in flux for clarity. 
The decomposition fitting is heuristic, and the model adopted for the fitting is composed of four 
Gaussian components, whereof two are intended for [O 1] AA6300, 6363 (shown as cyan lines), and 
the other one is intended for the component on the red side of [O I] lines which might originate 
from Ha, and/or [N 1n] \A6548, 6583 (shown as magenta lines), and the last component on the 
blue side and red side of the [O 1] lines (shown as the black lines) was added to improve the fitting. 


The two Gaussian components for [O 1] \A6300, 6363 share the same width and have flux ratio of 
f({OT]A6300) __ 
f(O 1|\6363) ~~ 
The fittings were performed for three spectra taken after the fast decline of the light curve, and one 


3. The fitting results shown as the red lines reproduce the observed spectra well. 


spectrum taken before the fast decline is shown for comparison. 
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Supplementary Information Figure 9: Spectral energy distribution (SED) of the new y-ray 
source. The energy range (100 MeV -300 GeV) is evenly divided in logarithmic space. The test 
statistic (TS) for each energy band is shown with the red bars, and the values are shown on the 
right axis. The blue points are the measured energy fluxes of sources with TS > 2 in individual 
energy bins, and the yellow points show the 95% confidence upper limits. The black solid line is 
the power law model with the dotted line as uncertainty. The top left panel is the SED for data in 
September and October 2022, where no source is detected in individual energy bands with signifi- 
cance above 95% confidence. The top right is the SED for data in November and December 2022, 
which shows the source is significantly detected and only detected in the 1—3 GeV energy band. 
The bottom panel is the SED for the data taken from September 1 to December 31. 
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Supplementary Information Figure 10: S-map grid of the field of SN 2022jli. Each stamp is 
built out of 2-month Fermi-LAT data of photons-yithin energy band 1 GeV to 3 GeV, with the year 
and months of the data shown in the left and top of the figure, e.g., the top right corner stamp 
corresponds to data taken between 2008 November | to 2009 January 1. The field of view of each 
stamp is 4° x 4°. The red plus symbol indicates the position of SN 2022jli, and the value in the 
bottom right corner of each stamp shows the highest S value with 0.25° radius of SN 2022jli. 
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Supplementary Information Figure 11: S-map grid of SN 2022)li field (continue of Supplemen- 


tary Information Fig|10). 75 
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Supplementary Information Figure 12: S-map grid of SN 2022jli field. The same as Supple- 
mentary Information Fig|10|except the energy of Photons is in the range of 300 MeV to 1 GeV. 
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Supplementary Information Figure 13: S map of SN 2022jli field (continue of Supplementary 


Information Fig|12) 77 
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Supplementary Information Figure 14: S-map grid of SN 2022jli field. The same as Supple- 
mentary Information Fig|10|except the energy of photons is in the range of 3 GeV to 10 GeV. 


March,April May,June  July,August Sep,Oct Nov,Dec 


2016 
a 
ew 


2017 


2021 2020 2019 2018 


2022 


2023 


20 
| 10 
0 


Supplementary Information Figure 15: S-map of SN 2022jli field (continue of Supplementary 


Information Fig|14) 79 
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Supplementary Information Figure 16: The Lomb-Scargle periodogram of y-ray photon arrival 
time. The Lomb-Scargle power spectra come from a series of unit values at the photon arrival time. 
The upper panel shows the result of photons in a narrow energy range of 1 GeV to 3 GeV, and the 
bottom panel shows the result from a broad energy range of 100 MeV to 300 GeV. The red line is 
for photons within 0.4° radius of the new y-ray source. The yellow lines are for photons within 
a specific radius of the 42 cataloged sources in 4FGL-DR3 within 10° radius of SN 2022jli. The 
radius for the 4FGL-DR3 source is chosen to have the same number of photons as the new y-ray 
source, i.e., 11 photons for 1 GeV to 3 GeV, and 37 photons for 100 MeV to 300 GeV. The black 
lines show the averaged power spectra of the 42 4FGL-DR3 sources, which provide information 
on the window function of the Fermi-LAT observation of the region. The period around 0.132 days 
(~ 190 minutes) in the averaged spectrum is clear, corresponding to typical observatory sky-survey 
profiles which run for two orbits, with the observatory rocked to the north for the first orbit and 
to the south for the second. As a result, the SN 2022jli region has been observed every two orbits 
of Fermi-LAT observation. There is no significant periodic signal at 12.4 days, which means the 
y-ray photon’s association to the 12.4-day period of SN 2022jli shown in Extended Data Fig. 
could not be due to the survey profile. The peaks around 12.4 days in the power spectra of the new 
y-ray source indicated by the vertical dashed line in red color, although not significant alone, are 
consistent with the result found in Extended Data Fig. 
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Supplementary Table 1. SN 2022jli photometry 


Epoch AB Magnitude Source Filter 
(JD-2,459,700) 


8.9154 15.366+0.053 ASAS-SN g 
12.6611 15.386+0.053 ASAS-SN g 
14.6775 15.602+0.171 ASAS-SN g 
16.6692 15.737+0.120 ASAS-SN g 
22.8739 15.477+0.060 ASAS-SN g 
25.6585 14.786+0.007 ATLAS o 
25.6657 14.775+0.007 ATLAS o 
25.6826 14.796+0.012 ATLAS o 
27.6503 16.1400.113 ASAS-SN g 
31.6665 16.245+0.078 ASAS-SN g 
33.0868 15.751+0.013 ATLAS c 
33.0953 15.696+0.011 ATLAS c 
33.1092 15.707+0.013 ATLAS c 
33.1134 15.718+0.017 ATLAS c 
33.6136 16.092+0.063 ASAS-SN g 
34.6683 16.324+0.125  ASAS-SN g 
35.6096 16.005+0.082 ASAS-SN g 
37.0850 15.568+0.011 ATLAS c 
37.1016 15.540+0.011 ATLAS c 
37.1042 15.920+0.072 ASAS-SN g 
37.1044 15.529+0.010 ATLAS c 
37.1150 15.579+0.018 ATLAS c 
37.6023 15.944+0.066 ASAS-SN g 
37.6553 14.867+0.007 ATLAS o 
37.6739 14.861 +0.008 ATLAS o 
37.6854 14.882+0.015 ATLAS o 
38.8515 15.774+0.051 ASAS-SN g 
41.5805 15.720+0.053 ASAS-SN g 
42.9147 15.482+0.039 ASAS-SN g 
45.0839 14.578+0.009 ATLAS o 
45.0959 14.517+0.009 ATLAS o 
45.1065 14.563 +0.008 ATLAS o 
45.1101 14.538 0.008 ATLAS o 
47.0924 14.489+0.009 ATLAS o 
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Supplementary Table 1 (contd) 


Epoch AB Magnitude Source Filter 
(JD-2,459,700) 


47.0965 14.520+0.009 ATLAS o 
47.1034 14.494+0.008 ATLAS o 
47.1159 14.496+0.019 ATLAS o 
47.8300 15.2750.047 ASAS-SN g 
49.8793 15.202+0.049  ASAS-SN g 
54.0836 14.482+0.005 ATLAS o 
54.0868 14.464+0.006 ATLAS o 
54.0966 14.463+0.005 ATLAS o 
54.1023 14.454+0.005 ATLAS o 
55.0878 14.434+0.006 ATLAS o 
55.0925 14.483 0.006 ATLAS o 
55.1004 14.475+0.006 ATLAS o 
55.1018 14.464+0.006 ATLAS o 
55.1105 14.455+0.006 ATLAS o 
55.1119 14.427+0.006 ATLAS o 
55.1183 14.48 10.009 ATLAS o 
55.1197 14.467+0.011 ATLAS o 
57.9530 15.244+0.024 ZTF g 
57.9530 15.244+0.024 ZTF g 
57.9540 15.220+0.029 ZTF g 
57.9540 15.220+0.029 ZTF g 
57.9540 15.220+0.029 ZTF g 
58.9335 15.217+0.035 ASAS-SN g 


Note. — Only part of the photometry is shown here for 


guidance regarding the form and content of the table. 
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Supplementary Table 2. Summary of Spectroscopic Observations of SN 2022jli 


Date (UTC) Phase (d) Telescope Spectrograph Exp (s) 


2022-06-29.423 55.254 P60 SEDM 1800 
2022-07-07.424 63.254 P60 SEDM 1800 
2022-07-19.453 75.284 P60 SEDM 1800 
2022-07-26.338 82.169 P60 SEDM 1800 
2022-08-04.418 91.249 P60 SEDM 1800 
2022-08-06.196 93.027 NOT ALFOSC 600 

2022-08-11.310 98.141 P60 SEDM 1800 
2022-08-16.288 103.119 P60 SEDM 1800 
2022-08-18.325 105.155 P60 SEDM 1800 
2022-08-20.0 106.830 P200 DBSP 300 

2022-08-23.371 110.202 P60 SEDM 1800 
2022-08-26.824 113.154 Magellan FIRE 507 

2022-08-28.254 115.084 P60 SEDM 1800 
2022-09-08.387 126.217 P60 SEDM 1800 
2022-09-16.216 134.046 P60 SEDM 1800 
2022-09-21.221 139.052 NOT ALFOSC 1200 
2022-10-02.305 150.136 P60 SEDM 1800 
2022-10-05.149 152.980 P60 SEDM 1800 
2022-10-11.159 158.989 P60 SEDM 1800 
2022-10-11.274 159.104 P200 DBSP 600 

2022-10-14.199 162.029 P60 SEDM 1800 
2022-10-21.179 169.009 P60 SEDM 1800 
2022-10-22.372 170.202 P60 SEDM 1800 
2022-10-30.337 178.168 P60 SEDM 1800 
2022-10-31.997 179.828 NOT ALFOSC 1200 
2022-11-02.0 180.830 P200 DBSP 420 

2022-11-15.133 193.964 P60 SEDM 1800 
2022-11-25.1 203.831 NOT ALFOSC 1200 
2022-11-27.78 205.909 P60 SEDM 1800 
2022-12-11.119 219.949 P60 SEDM 1800 
2022-12-15.137 223.967 Magellan IMACS 900 

2022-12-15.560 223.890 Magellan FIRE 2832 
2022-12-20.98 228.929 P60 SEDM 1800 


2022-12-27.91 235.921 MMT BINOSPEC 2400 
2023-01-12.77 251.907 MMT BINOSPEC 450 


83 


Supplementary Table 2 (cont’d) 


Date (UTC) Phase (d) Telescope Spectrograph Exp (s) 


2023-01-13.71 252.902 MMT BINOSPEC 1350 


2023-01-13.172 253.003 P60 SEDM 1800 
2023-01-14.50 253.881 VLT XSHOOTER 900 
2023-01-18.44 257.874 VLT XSHOOTER 900 
2023-01-21.45 260.875 VLT XSHOOTER 900 


2023-01-21.111 260.941 MMT BINOSPEC 1800 
2023-01-22.90 261.921 MMT BINOSPEC 2250 
2023-01-31.89 270.920 MMT BINOSPEC 1350 


2023-02-02.108 272.938 P60 SEDM 1800 

2023-02-03.829 274.660 NOT ALFOSC 1200 

2023-02-08.824 279.654 NOT ALFOSC 1200 
Note. — The phase is calculated with respect to the discovery date of 


JD= 2459704.67. 
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